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Abstract 

High-temperature series are computed for a generalized 3d Ising model with 
arbitrary potential. Two specific "improved" potentials (suppressing leading 
scaling corrections) are selected by Monte Carlo computation. Critical ex- 
ponents are extracted from high-temperature series specialized to improved 
potentials, achieving high accuracy; our best estimates are: 7 = 1.2371(4), 
v = 0.63002(23), a = 0.1099(7), 77 = 0.0364(4), (3 = 0.32648(18). By the 
same technique, the coefficients of the small-field expansion for the effective 
potential (Helmholtz free energy) are computed. These results are applied 
to the construction of parametric representations of the critical equation of 
state. A systematic approximation scheme, based on a global stationarity 
condition, is introduced (the lowest-order approximation reproduces the lin- 
ear parametric model). This scheme is used for an accurate determination 
of universal ratios of amplitudes. A comparison with other theoretical and 
experimental determinations of universal quantities is presented. 

Keywords: Critical Phenomena, Ising Model, High- Temperature Expansion, 
Critical Exponents, Critical Equation of State, Universal Ratios of Ampli- 
tudes, Effective Potential. 

PACS Numbers: 05.70.Jk, 64.60.Fr, 64.10.+h, 75.10.Hk, ll.10.Kk, 11.15.Me 



1 



Typeset using REVTeX 



I. INTRODUCTION 



According to the universality hypothesis, critical phenomena can be described by quan- 
tities that do not depend on the microscopic details of a system, but only on global prop- 
erties such as the dimensionality and the symmetry of the order parameter. Many three- 
dimensional systems characterized by short-range interactions and a scalar order parameter 
(such as density or uniaxial magnetization) belong to the Ising universality class. This 
implies that the critical exponents, as well as other universal quantities, are the same for 
all these models. Their precise determination is therefore important in order to test the 
universality hypothesis. 

The high-temperature (HT) expansion is one of the most effective approaches to the study 
of critical phenomena. Much work (even recently) has been devoted to the computation of 
HT series, especially for iV-vector models and in particular the Ising model. An important 
issue in the analysis of the HT series is related to the presence of non-analytic corrections 
to the leading power-law behavior. For instance, according to the renormalization group 
theory (see e.g. Ref. H), the magnetic susceptibility should behave as 



x = Cr 1 (1 + a ,it + a 0|2 r + ... + a hl t* + a h2 t 2 * + ... + a 2 ,i^ 2 + ...J , (1.1) 

where t = (T — T c )/T c is the reduced temperature. The leading critical exponent 7, and 
the correction exponents A, A 2 , are universal, while the amplitudes C and a^j are non- 
universal and should depend smoothly on any subsidiary parameter that may change T c , but 
does not affect the nature of the transition. Non-analytic correction terms to the leading 
power-law behavior, represented by noninteger powers of t, are related to the presence of 
irrelevant operators. For three-dimensional Ising-like systems the existence of leading cor- 
rections described by A ~ 0.5 is well established. In order to obtain precise estimates of the 
critical parameters, the approximants of the HT series should properly allow for the conflu- 
ent non-analytic corrections [0-BJ. The so-called integral approximants || can, in principle, 
allow for them (see e.g. Ref. [Kj for a review). However, they require long series to detect 



non-leading effects, and in practice they need to be biased to work well. Analyses meant to 
effectively allow for confluent corrections are generally based on biased approximants where 
the value of f3 c and the first non-analytic exponent A is given (see e.g. Refs. JTT|-|ToT|). It is in- 
deed expected that the leading non-analytic correction is the dominant source of systematic 
error. 

An alternative approach to this problem is the construction of a HT expansion where 
the dominant confluent correction is suppressed. If the leading non-analytic terms are not 
anymore present in the expansion, the analysis technique based on integral approximants 
should become much more effective, since the main source of systematic error has been elim- 
inated. In order to obtain an improved high-temperature (IHT) expansion, we may consider 
improved Hamiltonians characterized by a vanishing coupling with the irrelevant operator 
responsible for the leading scaling corrections. This idea has been pursued by Chen, Fisher, 
and Nickel ||, who studied classes of two-parameter models (such as the bcc scalar double- 
Gaussian and Klauder models). Such models interpolate between the spin-1/2 Ising model 
and the Gaussian model, and they are all expected to belong to the Ising universality class. 
The authors of Ref. || showed that improved models with suppressed leading corrections 



2 



to scaling can be obtained by tuning the parameters (see also Refs. [|3l,|T6f ) . This approach 
has been recently considered in the context of Monte Carlo simulations |T7|-pO|, using lat- 
tice 4 models. It is also worth mentioning the recent work [[H]] of Belohorec and Nickel 
in the context of dilute polymers, where a substantial improvement in the determination of 
the critical exponents v and uj was achieved by simulating the two-parameter Domb- Joyce 
model. 

We consider the class of scalar models defined on the simple cubic lattice by the Hamil- 
tonian 

n = -p £ + (1.2) 

<i,j> i 

where (3 = 1/T, < i,j > indicates nearest-neighbor sites, 0, are real variables, and V(4> 2 ) 
is a generic potential (satisfying appropriate stability constraints). The critical limit of 
these models is expected to belong to the Ising universality class (apart from special cases 
corresponding to multicritical points). Using the linked cluster expansion technique, we 
calculated the high-temperature expansion to 20th order for an arbitrary potential, general- 
izing the existing expansions for the standard Ising model (see e.g. Ref. |H| for a review of 
the existing HT calculations). In this work we will essentially consider and present results 
for a potential of the form 

U(0 2 ) = 2 + A 4 (0 2 - l) 2 + A 6 (0 2 - l) 3 ; (1.3) 

such a potential will be assumed in the following, unless otherwise stated. Within this 
family of potentials, improved Hamiltonians can be obtained by looking for values of the 
parameters A4 and Xq for which leading scaling corrections are suppressed. In particular we 
may keep A 6 fixed and look for the corresponding value A4 of A 4 that produces an improved 
Hamiltonian. Notice that, for generic choices of the Hamiltonian, A4 may not exist. This 
is the case of the O(N) 4 theory with nearest-neighbor couplings on a cubic lattice in the 
large- N limit, where it is impossible to find a positive value of A4 achieving the suppression 
of the dominant scaling corrections. Using the estimates of the leading scaling correction 
amplitudes reported in Ref. JT5], one can argue that the same is true for finite N > 3. As 



shown numerically by Monte Carlo simulations |17|- p0|1 , A4 exists in the case N = 1, which 
is the single-component <p A model (i.e. the model presented above with A6 = 0). By using 
finite-size techniques, Hasenbusch obtained a precise estimate of X\\ X\ = 1.10(2) |p0|| . m 
our work we will also consider the spin-1 (or Blume-Capel) Hamiltonian 

H = -P £ s iSj + Dj2sl (1.4) 

<i,j> i 

where the variables Sj take the values 0, ±1. In this case the value of D for which the leading 
scaling corrections are suppressed is D* = 0.641(8) p2| . 



From the point of view of the HT expansion technique, the main problem is the deter- 
mination of the improved Hamiltonian. Once the improved Hamiltonian is available, the 
analysis of its HT series leads, as we shall see, to much cleaner and therefore reliable results. 
A precise estimate of the parameters associated with an improved Hamiltonian is crucial in 
order to obtain a substantial improvement of the IHT results. As shown in Refs. P"7|.pU 
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Monte Carlo simulations using finite-size scaling techniques seem to provide the most effi- 
cient tool for this purpose. For comparison, in the case of the pure 4 theory (Ag = 0), our 
best estimate of A4 from the HT expansion is consistent with the above-mentioned Monte 
Carlo result A4 = 1.10(2), but it is affected by an uncertainty of about 10% (see Sec. |T|). 
So we decided to follow the strategy of determining the improved Hamiltonian by Monte 
Carlo simulations employing finite-size scaling techniques. For the 4 model, this work has 
been satisfactorily done by Hasenbusch [M. For the 6 model with Ae = 1, we performed 
Monte Carlo simulations in order to calculate A4, obtaining A4 = 1.90(4). 

The comparison of the results obtained from the three improved Hamiltonians considered 
strongly supports our working hypothesis of the reduction of systematic errors in the IHT 
estimates, and provides an estimate of the residual errors due to the subleading confluent 
corrections to scaling. 

The analysis of our 20th order IHT series allows us to obtain very precise estimates of the 
critical exponents 7, v and r\. Our estimates substantially improve previous determinations 
by HT or other methods. 

We extended our study to the small-field expansion of the effective potential, which is 
the Helmholtz free energy of the model. This expansion can be parametrized in terms of 
the zero-momentum n-point couplings g n in the symmetric phase. The analysis of the IHT 
series provides new results for the couplings g n , and leads to interesting comparisons with the 
estimates from other approaches based on field theory and lattice techniques. Moreover we 
improved the knowledge of the universal critical low-momentum behavior of the two-point 
function of the order parameter, which is relevant for critical scattering phenomena. 

By exploiting the known analytic properties of the critical equation of state, one may 
reconstruct the full critical equation of state from the small-field expansion of the effective 
potential, which is related to the behavior of the equation of state for small magnetization in 
the symmetric phase. This can be achieved by using parametric representations implement- 
ing in a rather simple way the known analytic properties of the equation of state. Effective 
parametric representations can be obtained by parametrizing the magnetization M and the 
reduced temperature t in terms of two variables R and 9, setting M oc R@6, t = R(l — 6 2 ), 
and H oc R^ & h{6). In this framework, following Guida and Zinn- Justin [P^j , one may develop 
an approximation scheme based on truncations of the Taylor expansion of the function h{9) 
around 9 = 0. Knowing a given number of terms in the small-field expansion of the effective 
potential, one can derive the same number of terms in the small-# expansion of h(6), with 
a dependence on an arbitrary normalization parameter p. One can try to fix p so that this 
small-6 1 expansion has the fastest possible convergence. We propose a prescription based on 
the global stationarity of the truncated equation of state with respect to the arbitrary pa- 
rameter p. This extends the stationarity condition of the linear model (i.e. the lowest-order 
non-trivial approximation) discussed in Refs. [24-27fl. Using the IHT results for 7, v and the 
first few coefficients of the small-field expansion of the effective potential, we constructed 
approximate representations of the full critical equation of state. From them we obtained 
accurate estimates of many ratios of universal amplitudes. Varying the truncation order of 
h{9), we observed a fast convergence, supporting our arguments. 

For our readers' convenience, we collected in Table [Xllfj a summary of all the results 
obtained in this paper. There they can find new estimates of most of the universal quan- 
tities (exponents and ratios of amplitudes) introduced in the literature to describe critical 
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phenomena in 3c? Ising-like systems. The only important quantity for which we have not 
been able to give a good estimate is the exponent u, which is related to the leading scaling 
corrections. We mention that a precise estimate of u has been reported recently in Ref. 
20fl : uj = 0.845(10). It has been obtained by a Monte Carlo study using a finite-size scaling 



method. 

The paper is organized as follows. 

In Sec. [TT] we discuss the main features of improved Hamiltonians from the point of view 
of the renormalization group. 

In Sec. |TT| we describe our Monte Carlo simulations and present the estimates of A4 for 
the potential ( |1.3| ) with A 6 = 1. 

Sec. |V| is dedicated to the determination of the critical exponents from the IHT expan- 
sion. We present estimates of all relevant critical exponents (except for uj), and compare 
our results with other theoretical approaches and experiments. 

In Sec. [V] we study the small-field expansion of the effective potential. We present 
estimates of the first few coefficients of the expansion. We discuss the relevance of the 
determination of the zero-momentum four-point renormalized coupling for field-theoretical 
approaches (Sec. |VB|) . 

Sec. |VT] presents a study of the low-momentum behavior of the two-point function in 
the critical region. Estimates of the first few coefficients of its universal low-momentum 
expansion are given. 

In Sec. |VII| we study the critical equation of state, which gives a description of the 
whole critical region, including the low-temperature phase. Using the estimates of the 
critical exponents and of the first few coefficients of the small-field expansion of the effective 
potential, the critical equation of state is reconstructed employing approximate parametric 
representations (Sec. |V 11 A| ). In Sec. |V11B| we present our method, based on the global 
stationarity of the approximate equation of state. Relevance to the e-expansion is discussed 
in Sec. |V 11 C[ In Sec. |VllD| we apply the results of Sec. |VllB| to the computation of universal 
ratios of amplitudes, using as inputs the results of the IHT expansion. The results are then 
compared with other theoretical estimates and with experimental determinations. For sake 
of comparison, we also present results for the two-dimensional Ising model. 

Many details of our calculations are reported in the Appendices. App. [A] contains a 
detailed description of our HT calculations, i.e. the list of the quantities we have considered 
and the description of the method we used to generate and analyze the HT series. We 
report many details and intermediate results so that the reader can judge the quality of 
the results we will present. In App. || we present the notations for the critical amplitudes, 
and report the expressions of the universal ratios of amplitudes in terms of the parametric 
representation of the critical equation of state. In App. |C| we discuss in more detail the 
approximation scheme for the parametric representation of the equation of state based on 
stationarity. 



II. IMPROVED HAMILTONIANS 

As discussed in the introduction, we will work with "improved" Hamiltonians, i.e. with 
models in which the leading correction to scaling has a vanishing (in practice very small) 
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amplitude. 

To clarify the basic idea, let us consider a model with two relevant operators (the thermal 
and the magnetic ones) and one irrelevant operator. If r, k, and fi are the associated non- 



linear scaling fields, the singular part of the free energy F sing has the scaling form [28[] 



F smg (r, K ,/i) = \r\ dv f ± {K\r\-^ 2 -^ 2 ^\rf) . (2.1) 

where the function f± depends on the phase of the model. Since the operator associated 
with /j, is irrelevant, A is positive and /i|r| A — > at the critical point. Therefore we can 
expand the free energy obtaining 

oo 

F smg (r,K,/i) = |r| d ^/ n , ± ( K |rr( rf+2 -^/ 2 )^|rr A . (2.2) 



n=0 



The presence of the irrelevant operator induces non-analytic corrections proportional to 
|r| nA . Now, let us suppose that the Hamiltonian of our model depends on three parameters 
r, h, and A, where r is associated to the temperature, h is the magnetic field and A is an 
irrelevant parameter. For each value of A and for h = 0, the theory has a critical point for 
r = r c (A). The non-linear scaling fields r, k, and [i are analytic functions of the parameters 
appearing in the Hamiltonian, and therefore we can write 

t = t + t 2 g lT {\) + h 2 g 2r (X) + 0(t 3 , th 2 , /i 4 ), (2.3) 
K = h[l + tg lK {\) + h 2 g 2K (\) + 0(t 2 , th 2 , h 4 )] , (2.4) 
fi = gifx{\) + tg 2lx {\) + h 2 g^(X) + 0(t 2 , th 2 , h 4 ), (2.5) 

where t = r — r c (A). Substituting these expressions into Eq. (|2~2|), we see that, if gi^{\) ^ 0, 
the free energy has corrections of order t nA . For the susceptibility in zero magnetic field we 



obtain the explicit formula [29 



X=t-~< £ XlMW 1 ^ 71 +t'~ a E X2, m n(A)t mA+ ™ + E X3,n(A)t™, (2.6) 
m,n=0 m,n=0 n=0 

where the contribution proportional to t 1_Q stems from the terms of order h 2 appearing 
in the expansion of r and /i, and the last term is the contribution of the regular part of 
the free energy. Notice that it is often assumed that the regular part of the free energy 
does not depend on h. If this were the case, we would have X3,n(A) = 0. However, for the 
two-dimensional Ising model, one can prove rigorously that X3,o PD|J31|, showing the 



incorrectness of this conjecture. For a discussion, see Ref. p2 |. 

In many interesting instances, it is possible to cancel the leading correction due to the 
irrelevant operator by choosing A = A* such that gi^A*) = 0. In this case fir A ~ t 1+A , so 
that no term of the form t mA+n , with n < m, will be present. In particular the leading term 
proportional to t A will not appear in the expansion. 

In general other irrelevant operators will be present in the theory, and therefore we expect 
corrections proportional to t p with p = n\ + n 2 A + J2i m iAi, where Aj are the exponents 
associated to the additional irrelevant operators. For A = A* the expansion will contain only 
terms with rii >n 2 . 
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It is important to notice that, by working with A = A*, we use a Hamiltonian such that 
the non- linear scaling field /i vanishes at the critical point. This property is independent of 
the observable we are considering. Therefore all quantities will be improved, in the sense 
that the leading correction to scaling, proportional to t A , will vanish. We will call the 
Hamiltonians with A = A* "improved Hamiltonians" . 



III. DETERMINATION OF THE IMPROVED PARAMETERS 



The Hamiltonian defined by Eqs. ( |1.2| ) and (|1.3| ) with A6 = was considered in Ref. |20 
where it was shown that the leading correction to scaling cancels for A4 = 1.10(2). Here 
we will also consider the case Xq = 1, and determine the corresponding A4 using a method 



similar to the one discussed in Ref. |L7H 



The idea is the following: consider a renormalizat ion-group invariant observable O on a 
finite lattice L and let O* be its value at the critical point, i.e. 

O* = lim lim W,X A ,L). (3.1) 

The quantity O* is a universal number and therefore it will be independent of A4. The 
standard scaling arguments predict 

0(&(A 4 ), A 4 , L) « O* + ai {\ A )L-" + a 2 (\ 4 )L~ 2 " + ... + 6 1 (A 4 )L- 2 . . . (3.2) 

where u = A/V, u 2 = A 2 /z/, A 2 being the next-to-leading correction-to-scaling exponent. 
Since for A 4 = A 4 , ai(A|) = e^A^) = . . . = 0, for A 4 ~ A4 we can rewrite the previous 
equation as 

0((3 C (X 4 ), A 4 , L) « O* + (A 4 - A!) (anL-* + a 21 L- 2 " + ...)+ b^XDL^ . . . (3.3) 

Now, suppose we know the exact value O*, and let us define A| ff (L) as the solution of the 
equation 

0((3 c (\f(L)),\f(L),L) = 0*. (3.4) 
From Eq. ( |3.3|) we obtain immediately 

Af (L) = X* 4 - ^r- W2 + . . . (3.5) 

Since u 2 > lu, X'f(L) converges to A4 as L — » 00. For the 3d Ising universality class, oj 2 — 2u 
H,H and to ~ 0.85 |C 



In order to apply this method in practice we need two ingredients: a precise determination 
of /3 C (A 4 ) and an estimate of O* . 

Very precise estimates of /3 C (A 4 ) can be obtained from the analysis of the HT series of 
the susceptibility Xi that we have calculated to 0(f3 20 ). For A6 = and 1.0 < A 4 < 1.2 the 
values of /3 C (A 4 ) can be interpolated by the polynomial 

&(A 4 ) = 0.40562043 + 0.00819000 A 4 - 0.04626355 X\ + 0.01235674 A^ ± 0.0000014. (3.6) 
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In particular, for A4 = 1.10, we have /5 C (1.10) = 0.3750973(14), to be compared with 



&(1.10) = 0.3750966(4) of Ref. ||. For A 6 = 1 and 1.8 < A 4 < 2.0 — as we shall see 



this is the relevant interval — we have 

/3 C (A 4 ) = 0.68612192 - 0.18274273 A 4 + 0.02634688 \\ - 0.00102710 A^ ± 0.0000018. (3.7) 

The second quantity we need is an observable O such that O* can be computed with 
high precision. We have chosen the Binder parameter 



(m 4 ) ( 
Q = 1 — ^2,1 (3-8) 



where m is the magnetization. A precise estimate of Q was obtained in Ref. ]17J by means 
of a large-scale simulation of the spin-1 model. They report 

Q* = 0.62393(13+35+5), (3.9) 

where the error is given as the sum of three contributions: the first is the statistical error; 
the second and the third account for corrections to scaling. We have tried to improve this 
estimate by performing a high-precision Monte Carlo simulation^ of the Hamiltonian ( |1.2|) 
for A 6 = and by computing Q for A 4 = 1.10, which is the best estimate of A4. Since this 
Hamiltonian is improved (i.e. the leading correction to scaling vanishes), we expect to be 
able to obtain a reliable estimate from simulations on small lattices for which it is possible to 
accumulate a large statistics. The results^ are reported in Table |. There are three different 
sources of error: we report the statistical error, the variation of the estimate of the Binder 
parameter when A 4 varies within the interval 1.08-1.12 (due to corrections to scaling of order 
L~^, which are not completely suppressed, since the value used for A 4 is not exactly equal to 
A 4 ), and the variation of Q when f3 c varies within one error bar. The values of L we use are 
relatively small (L < 12) and one could be afraid that next-to-leading corrections still give 
a non-negligible systematic deviation. Our data do not show any evidence of such an effect, 
and the estimates for different values of L are consistent. Using the estimate obtained for 
L = 12, we get the final result 

Q* = 0.62388(32) (3.10) 

(the uncertainty is obtained assuming independence of systematic and statistical errors), 
which is in agreement with the estimate (|3.9|) with a slightly smaller error bar. 



l We used the Brower-Tamayo algorithm |35|| , each iteration consisting of a Swendsen-Wang update 
of the sign of <j) and of a Metropolis sweep. 

2 The simulations were performed before the appearance of Ref. |2(| and the generation of the HT 
series, when only a very approximate expression for A4 existed. Therefore the runs were not made 
at the correct values of A 4 and (3 C . The values reported in Table | have been obtained from the 
Monte Carlo data by means of a standard reweighting technique. 
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TABLE I. For several values of lattice size L and for \q = 0, we report: the values of the 
parameters used in the Monte Carlo simulation A4 irun , /3 run , the number of Monte Carlo iterations 
./Vi ter , each iteration consisting of a standard Swendsen-Wang update and of a Metropolis sweep, 
and the the estimate of the Binder parameter Q at A = 1.10, (3 = 0.3750973. The reported error 
on Q is the sum of three terms: the statistical error, and the errors due to the uncertainty of A* 
and /3 C (A). 

L A4 jrU n ftun -Alitor Q 

~6 1.100 0.375 91 x 10 6 0.62370(15+37+2) 

7 1.100 0.375 78 x 10 6 0.62386(17+33+2) 

9 1.080 0.376 178 x 10 6 0.62389(12+26+3) 

12 1.105 0.375 305 x 10 6 0.62387(10+25+5) 



TABLE II. For several values of lattice size L and for A6 = 1, we report: the values of the 
parameters used in the Monte Carlo simulation A4 jrun , /3 run , the number of Monte Carlo iterations 
iViten each iteration consisting of a standard Swendsen-Wang update and of a Metropolis sweep, 
and the estimate of A eff (L), the solution of Eq. ( |3.4| ) . The error is reported as the sum of three 
terms: the statistical error on Q, the error on Q* , and the uncertainty of the critical value (3 C {X). 



L 


A4,run 


/^run 


Mter 


Af(£) 


6 


1.900 


0.427 


100 x 10 6 


1.915(9+19+1) 


7 


1.900 


0.427 


252 x 10 6 


1.894(6+21+2) 


9 


1.920 


0.425 


214 x 10 6 


1.897(9+25+3) 


12 


1.920 


0.425 


294 x 10 6 


1.904(9+32+6) 



We have next determined A4 for the model with Hamiltonian ( IOJ ) and Xq = 1 using 
the method presented above. Estimates of Xf^(L) are reported in Table [D], from which we 
conclude 

A* = 1.90(4). (3.11) 

Notice that the last three points show a small upward trend which, although consistent with 
a statistical effect, could be a systematic increase^ due to the corrections of order £-< 4, 2+ £ + 
With the present statistical errors we cannot distinguish between these two possibilities: we 
have considered as our final estimate the average of the results obtained for L = 7,9, 12, 
but we cannot exclude that the correct A4 is slightly larger than our estimate. However the 
quoted error should be large enough to include this systematic increase. 



3 This effect could also be due to the fact that Q* is only approximately known. To check if this is 
the case, we have computed \f [ ±(L), the solutions of Eq. (3.4) with the r.h.s. replaced by Q* ±ctq, 
gq being the error on Q* . If the increase is associated with the uncertainty of Q* , we should 
observe that A| f ! t (L) have opposite trends, one increasing, the other decreasing. In the present 
case A| f j_(L) are both increasing, thereby excluding that this effect is due to the uncertainty of Q*. 
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As explained in the introduction, the analysis of HT series for the determination of 
universal quantities is sensitive to non-analytic scaling corrections. As we will discuss below, 
one can use this fact to obtain a rough estimate of the optimal value of A. 

Consider for example the zero-momentum four-point coupling constant #4 defined by 

^ = (3-12) 

where Xi £ 2 an d X4 are > respectively, the magnetic susceptibility, the second-moment corre- 
lation length and the zero-momentum four-point connected correlation function (definitions 
can be found in App. [A 1|) . We have chosen this observable because it appears to be affected 
by large corrections to scaling, but the method can be applied to any universal quantity. 
From the discussion of Sec. [T| we have for (3 — > (3 C 

gM = 9l + CA(Pc-P) A + ... (3.13) 

where g\ is a universal constant, and ca is a non- universal amplitude depending on the 
Hamiltonian. For improved models, as discussed before, ca = 0. The traditional methods 
of analysis, e.g. those based on Pade (PA) and Dlog-Pade (DPA) approximants, are unable 
to handle an asymptotic behavior like (|3.13|) unless A is an integer number, thus leading to 
a systematic error. Integral approximants allow for non-analytic scaling corrections, but, as 
already said, with the series of moderate length available today, they need to be biased to 
give correct results: without any bias they give estimates that are similar to those of PA's 



and DPA's |14|]. At present the only analyses that are able to effectively take into account 
the confluent corrections use biased approximants, fixing the value of j3 c and of the first 
non-analytic exponent A (see e.g. Refs. [p!0|- [T^ , p^JT5| , |36| for a discussion of this issue and 



for a presentation of the different methods used in the literature). The method we use has 



been proposed in Ref. [11 1 and generalized in Ref. [12]. The idea is to perform a Roskies 



transform (RT), i.e. the change of variables 

z = l-(l-/?/&) A , (3.14) 

so that the non-analytic terms in (3 C — (3 become analytic in 1 — z. Therefore the analysis 
of the resulting series by means of standard approximants should give correct results. For 
the models we are considering the exponent A is approximately 1/2 (e.g., Ref. |20|] reports 
A = 0.532(6)); for simplicity we have used the transformation ( |3.14| ) with A = 1/2. 

We have analyzed the HT expansion of g± for the model with Hamiltonian (|1.2|) and 
A 6 = for several values of A 4 . We computed PA's, DPA's, and first-order integral (IAl) 
approximants of the series in (3 and of its Roskies transform (RT) in z. In Fig. [I] we plot 
the results as a function of Aj . The reported errors are related to the spread of the results 
obtained from the different approximants, see App. | A 4| for details. The estimates obtained 



from the RT'ed series are independent of A4 within error bars, giving the estimate g\ ~ 23.5, 
in agreement with previous analyses of the HT expansion of the Ising model on various 



lattices using the RT or other types of biased approximants |1^JT5|J36]| (in Sec. [VB| we will 
improve this estimate by analyzing the IHT expansion). The independence of the result 
from the value of A4 clearly indicates that the RT is effectively able to take into account the 
non-analytic behavior ( |3.13| ). On the other hand the analysis of the series in (3 gives results 
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FIG. 1. Comparison of the determination of 54 (plotted vs. I/A4) from HT series without 
(direct) and with the Roskies transform (RT) for the pure (/> 4 lattice model. The dashed line marks 
the more precise estimate (with its error) we derived from the analysis of the IHT expansion. 



which vary with A4 more than the spread of the approximants: for instance, the analysis of 
the series of the standard Ising model, corresponding to A4 = 00, gives results that differ 
by more than 5% from the estimate quoted above, while the spread of the approximants is 
much smaller. Clearly there is a large systematic error. It is important to notice that the 
direct analysis and the RT one coincide when 1.0 ?S A4 ^$ 1.2, i.e. in the region in which the 
leading non-analytic corrections are small. This fact confirms our claim that the observed 
discrepancies are an effect of the confluent corrections. 

The results presented above can be used to obtain an estimate of A4 from the HT series 
alone: A4 should fall in the interval in which the direct analysis gives results compatible 
with those obtained from the RT'ed series. As we already mentioned, for A6 = we obtain 
A4 = 1.1(1), while for \ e = 1 we get A4 = 1.9(1). The latter estimate was indeed the starting 
point of our Monte Carlo simulation. We have also tried to estimate A4 in more direct ways, 
but all methods we tried were even less precise. 

A similar method for the determination of the improved Hamiltonian from the HT series 
was presented in Ref. ||. The optimal value of the parameter (called y in Ref. ||) was 
determined comparing the results for the critical point j3 c {y) obtained using IAl's and DPA's: 
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y* is estimated from the value at which DPA and IA1 estimates of f3 c (y) agree between 
each other. It should be noticed that, for the double-Gaussian model, partial differential 
approximants and a later analysis of Nickel and Rehr |T6] using a different method gave 
significantly different estimates of y*. 



IV. CRITICAL EXPONENTS 



The analysis of Sec. |T| is encouraging and supports our basic assumption that the system- 
atic error due to confluent singularities is largely reduced when analyzing IHT expansions. 
To further check this hypothesis we will compare results obtained from different improved 
Hamiltonians. This will provide an estimate of the remaining systematic error which is not 
covered by the spread of the results from different approximants. 

The definition of the quantities we have considered and a detailed description of the 
method we used to generate and analyze the HT series is presented in App. [A]. 

We computed (5 C and 7 from the analysis of the HT expansion of the magnetic suscepti- 
bility. We cosidered integral approximants of first, second and third order. After a careful 
analysis we preferred the second-order integral approximants (IA2's), which turned out to 
be the most stable: most of the results we present in this section and in the related App. |A 3| 
have been obtained by using IA2's. As a further check of the effectiveness of the approxi- 
mants employed, we made use of the fact that \ (and £ 2 ) must present an antiferromagnetic 



singularity at [3f = —f3 c of the form [[!7|] 

X = c + c 1 (/3-/3 c af ) 1 " a + ..., (4.1) 

where a is the specific heat exponent, Cj are constants and the ellipses represent higher-order 
singular or analytic corrections. We verified the existence of a singularity at (3 ~ — f3 c in the 
approximants, and calculated the associated exponent. We also considered approximants 
that were biased requiring the presence of two symmetric singularities at (3 = ±/3 c f|; the 
results obtained are consistent with the predicted behavior ( f4~T| ) (see App. |A3| and related 
Tables) . 

The exponent v was obtained from the series of the second-moment correlation length 

e = ^-(P c -pr, (4.2) 

where are the moments of the two-point function. We followed the procedure suggested 
in Ref. [J38||, i.e. we used the estimate of (5 C obtained from x t° bias the analysis of £ 2 . For 



this purpose we used IA's biased to f3 c . We also considered approximants biased to have a 
pair of singularities at ±/3 c . 

In Table [III] we report the results obtained for the Hamiltonians ( |1.2j ) with \ = 0, 
A 4 = 1.10 and with A 6 = 1, A 4 = 1.90, and for the Hamitonian Q with D = 0.641. The 
errors are given as a sum of two terms: the first one is computed from the spread of the 
approximants; the second one is related to the uncertainty of the value of A4 and D*, and 
it is evaluated by changing A4 in the range 1.08-1.12 for X e = and 1.86-1.94 for X e = 1, 
and D in the range 0.633-0.649 for the spin-1 model. There is a good agreement among the 
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TABLE III. Our final estimates of 7, r/ and a. The error is reported as a sum of two 
terms: the first one is related to the spread of the approximants; the second one is related to the 
uncertainty of the value of A| . 





7 


V 


V 


a 


A 6 = 


1.23732(24+16) 


0.63015(13+12) 


0.0364(3+1) 


0.0213(13+1) 


A 6 = l 


1.23712(26+31) 


0.63003(13+23) 


0.0363(3+2) 


0.0213(14+2) 


spin-1 


1.23680(30+12) 


0.62990(15+8) 


0.0366(3+1) 


0.0202(10+1) 



estimates of 7 and v obtained from the three improved Hamiltonians considered. This is 
an important check of our working hypothesis, i.e. that systematic errors due to confluent 
corrections are largely reduced. This will be also confirmed by the results for the universal 
ratios of amplitudes. We determine our final estimates by combining the results of the three 
improved Hamiltonians: as estimate we take the weighted average of the three results, and 
as estimate of the uncertainty the smallest of the three errors. We obtain for 7 and v 

7 = 1.2371(4), (4.3) 
v = 0.63002(23), (4.4) 

and by the hyperscaling relation a = 2 — 3u 

a = 0.1099(7). (4.5) 



In App. |A 3| we also report some further checks using the Monte Carlo estimate of j3 c reported 
in Ref. |2tJ to bias the analysis of the series. The results are perfectly consistent. We mention 



that from the analysis of the antiferromagnetic singularity, cf. Eq. (|4. 1|) , we obtain the 
estimate a = 0.105(10), which is consistent with result Q4.5|) obtained assuming hyperscaling. 

From the results for 7 and u, we can obtain rj by the scaling relation 7 = (2 — r\)v. This 
gives 77 = 0.0364(10), where the error is estimated by considering the errors on 7 and v as 
independent, which is of course not true. We can obtain an estimate of 77 with a smaller, yet 
reliable, error using the so-called critical point renormalization method (CPRM) (see Ref. 
and references therein). We obtain the results reported in Table [TTI|, with considerably 
smaller errors. Our final estimate is 

77 = 0.0364(4). (4.6) 
Moreover using the scaling relations we obtain 

5 = = 4.7893(22), (4.7) 

(3 = V - (1 + v ) = 0.32648(18) (4.8) 

(the error on j3 has been estimated by considering the errors of v and 77 as independent). 

Finally we consider the universal critical exponent describing how the spatial anisotropy, 
which is present in physical systems with cubic symmetry (e.g. uniaxial magnets), vanishes 
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when approaching the rotationally-invariant fixed point |39[]. For this class of systems the 
two-point function G(x) is not rotationally invariant. Therefore non-spherical moments are 
in general non-vanishing, but near the critical point they are depressed with respect to 
spherical moments carrying the same naive physical dimensions by a factor £~ p , where p is a 
universal critical exponent. From a field-theoretical point of view, space anisotropy is due to 
non rotationally invariant irrelevant operators in the effective Hamiltonian, whose presence 
depends essentially on the symmetries of the physical system, or of the lattice formulation. 
In Table |TXT| we report the results for a = 2 — p as obtained by analyses of the first non- 
spherical moments (cf. Eq. ( |A2| )) using the CPRM. The exponent a turns out to be very 
small: 



a = 0.0208(12), 



(4.9) 



and p = 1.9792(12). 

In Table |V] we compare our results with some of the most recent estimates of the critical 
exponents 7, u, r], a, and (3. The table should give an overview of the state of the art for 
the various approaches. Let us first note the good agreement of our IHT estimates with 
the very precise results of the recent Monte Carlo simulations (MC) of Refs. [ |iT|pi?| , [2"0"| . The 
small difference with the HT estimates of Ref. |T3 (obtained from the standard Ising model) 
may be explained by the difficulty of controlling the effects of the confluent singularities, 
and by a systematic error induced by the uncertainty on the external input parameters ((3 C 
and A) that are used in their biased analysis. The estimates of Refs. |5|,|7|,§,|n| have been 
obtained from a HT analysis of two families of models, the Klauder and the double-Gaussian 
models on the bcc lattice. The results of these analyses are in good agreement with our IHT 
estimates, especially those by Nickel and Rehr pi. The HT series for the double-Gaussian 
model were analyzed also in Ref. || where a higher estimate of 7 was obtained. As pointed 
out in Ref. |T6|] , the discrepancy is essentially due to the use of a higher estimate of the 
improvement parameter y* with respect to that used in Ref. |16j (see the discussion at the 
end of Sec. p| ). Refs. |]8|JT6|1 report also estimates of a obtained analyzing the singularity of 
the susceptibility at the antiferromagnetic critical point. The result agrees with our estimate. 
The agreement with the field-theoretical calculations is overall good. The slightly larger 

(using 0(g 7 ) 



BO 



series p6|,|50|) may be 



result for 7 obtained in the analyses of Refs. 
due to an underestimate of the systematic error due to the non-analyticity of the Callan- 
Symanzik /5-function. Similar results have been obtained by Kleinert, who resummed the 



0(g ) expansion by a variational method jEJ, still neglecting confluent singularities at the 



infrared- stable fixed point. We shall return on this point later. A better agreement is found 
with the analysis of the d=3 (^-expansion performed by Murray and Nickel, who allow for 
a more general non-analytic behavior of the /3-function ptj||. In Table |V|, we quote two 



errors on the results of Ref. ||6| : the first one is the resummation error, and the second one 
takes into account the uncertainty of g*, which is estimated to be ~ 0.01. The results of 
the e-expansion were obtained from the 0(e 5 ) series calculated in Refs. [pI,S2]|. We report 
estimates obtained by performing standard analyses (denoted as "free") and constrained 
analyses (denoted by "be" ) that incorporate the knowledge of the exact two-dimensional 
values. Both are essentially consistent with our IHT estimates, but present a significantly 
larger uncertainty. In Table [TV] we also report the results obtained by approximately solving 
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TABLE IV. Theoretical estimates of critical exponents. See text for 
in the first column. For values marked with an asterisk, the error is not 
reference. 



explanation of symbols 
quoted explicitly in the 
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V 


V 


a (3 


IHT 


1.2371(4) 


0.63002(23) 


0.0364(4) 


0.1099(7) 0.32648(18) 



HT (sc) H 
HT (bcc 
HT " 
HT 
HT 
HT 
HT 
HT 



m 



MC H 



1.2388(10) 

1.2384(6) 

1.237(2) 

1.239(3) 

1.2395(4) 

1.2378(6) 

1.2385(15) 

1.2385(25) 



0.6315(8) 
0.6308(5) 
0.6300(15) 
632+ - 002 

U.OO^_ QQ3 

0.632(1) 
0.63115(30) 

0.6305(15) 



0.0359(7) 



11(2) 
105(7) 



mc yj 
mc g 



mc gg| 

mc m 



mc m 

mc g 



1.2367(11) 



1.237(2) 



0.6296(7) 
0.6298(5) 
0.6294(10) 
0.6308(10) 

0.625(1) 
0.6301(8) 



0.0358(9) 
0.0366(8) 
0.0374(12) 



0.025(6) 
0.037(3) 



0.110(2) 



0.3269(6) 
0.3267(10) 



e-exp. frec g 
e-exp. bc 
e-exp. bc [II 



d=3 g-exp. |4| 
d=3 g-exp. [45 
d=3 g-exp. [46 
d=3 g-exp. m 



1.2355(50) 
1.2380(50) 
1.240(5) 



0.6290(25) 
0.6305(25) 
0.631(3) 



0.0360(50) 
0.0365(50) 



0.3257(25) 
0.3265(15) 



1.2396(13) 
1.241* 

1.2378(6+18) 
1.2405(15) 



0.6304(13) 
0.6305* 
0.6301(5+11) 
0.6300(15) 



0.0335(25) 
0.0347(10) 
0.0355(9+6) 
0.032(3) 



0.109(4) 0.3258(14) 



ERG [48| 




0.618(14) 


0.054* 


ERG |4S] 


1.247* 


0.638* 


0.045* 



the exact renormalizat ion-group equation (ERG); they seem to be much less precise. A 
more complete list of references pertaining to the theoretical determination of the critical 



exponents can be found in Ref. j4"4* |. Concerning the exponent a related to the rotational 
symmetry, the IHT results represent a substantial improvement of the estimates obtained 
by various approaches (HT and field theory) presented in Ref. [[^J. 

Experimental results have been obtained studying the liquid-vapor transition in sim- 
ple fluids, and the different critical transitions in multicomponent fluid mixtures, uniaxial 
antiferromagnetic materials and micellar systems. Many recent estimates can be found in 
Refs. [43j54,|55|. In Table |V| we report some experimental results, most of them published 
after 1990. It is not a complete list of the published results, but it may be useful to get an 
overview of the experimental state of the art. Even if the systems studied are quite different, 
the results substantially agree, although, looking in more detail, as already observed in Ref. 
one can find small discrepancies. Moreover they substantially agree with the theoretical 
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TABLE V. Experimental estimates of critical exponents, lv denotes the liquid-vapor transition 
in simple fluids, bm refers to a binary fluid mixture, ms to a uniaxial magnetic system, and mi to 
a micellar system. 



Ref. 7 v rj a f3 



lv 


a o.no5±™ 

57| 0.1075(54) 

H 0.1084(23) 

59| 0.111(1) 0.324(2) 

6(| 0.341(2) 

61| 0.042(6) 


bm 


6| 0.104(11) 

6| 1.233(10) 0.327(2) 

6| 1.09(3) 

6| 1.26(5) 0.64(2) 

66j 1.24(1) 0.606(18) 0.077(44) 0.319(14) 
67| 0.105(8) 

6| 0.324(5), 0.329(2) 
6§ 0.329(4), 0.333(2) 
7TJ 0.610(6) 

7| 0.336(30) 
7| 1.228(39) 0.628(8) 0.0300(15) 


ms 


7| 1.25(2) 0.64(1) 

7| 0.115(4) 0.331(6) 

7§ 0.11(3) 

76j 0.325(2) 

7^ 0.11(3) 

7| 1.25(2) 0.315(15) 


mi 


7| 0.34(8) 

H 1.18(3) 0.60(2) 

8| 1.216(13) 0.623(13) 0.039(4) 

H 1.237(7) 0.630(12) 

83| 1.25(2) 0.63(1) 

H 1.17(11) 0.65(4) 



predictions discussed above, confirming the fact that all these transitions are in the Ising 
universality class. It should also be noticed that the experimental results are less accurate 
than the theoretical estimates. 

V. THE EFFECTIVE POTENTIAL 
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A. Small-field expansion of the effective potential in the high-temperature phase 



The effective potential (Helmholtz free energy) is related to the (Gibbs) free energy of 
the model. Indeed, if M = ((f)) is the magnetization and H the magnetic field, one defines 

F{M) = MH — ^- log Z(H), (5.1) 

where Z(H) is the partition function and the dependence on the temperature is always 
understood in the notation. 

The global minimum of the effective potential determines the value of the order parame- 
ter which characterizes the phase of the model. In the high-temperature or symmetric phase 
the minimum is unique with M — 0. According to the Ginzburg-Landau theory, as the tem- 
perature decreases below the critical value, the effective potential takes a double-well shape. 
The order parameter does not vanish anymore and the system is in the low-temperature or 
broken phase. Actually in the broken phase the double-well shape is not correct because the 
effective potential must be convex f85fl . In this phase it should present a flat region around 



the origin. 

In the high-temperature phase the effective potential admits an expansion around M = 0: 

OO 1 

A jr = HM) _ m = £ __„. i; ,U^ (5.2) 

The coefficients a 2 j can be expressed in terms of renormalization-group invariant quantities. 
Introducing a renormalized magnetization 

2 t(t,H = oyM(t,H? 
V = *t,H = 0) • (5 ' 3) 

where t is the reduced temperature, one may write 

AT = ImV + X'"'' r " ^^r 1 '. (5.4) 

Here m = l/£, g 2 j are functions of t only, and d is the space dimension. In field theory ip 
is the expectation value of the zero-momentum renormalized field. For t — > the quantities 
92j approach universal constants (which we indicate with the same symbol) that represent 
the zero-momentum 2j-point renormalized coupling constants. By performing a further 
rescaling 

m (d~2)/2 

ip = z (5.5) 

\[9i 



in Eq. Q5.4Q , the free energy can be written as 



m d 



AT= — A(z), (5.6) 

94 
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where 



and 



r 2j — 3 > 3- (5.8) 

04 

One can show that z oc t~*M, and that the equation of state can be written in the form 

H oc t^4^. (5.9) 

az 

The effective potential JF(M) admits a power-series expansion also near the coexistence 
curve, i.e. for t < and H = 0. If M = lim iJ _ >0 + M(H), for M > M (i.e. for # > 0) we 
have 

5.F = JF(M) - JF(M ) = \ a A M ~ M oY- (5-10) 
In terms of the renormalized magnetization if we can rewrite 

sr = \m\f - v, ) 2 + J2^ 3{d ~ 2)/2 \gjbp - VoY, (5.11) 

where m = and £ _ is the second-moment correlation length defined in the low- 

temperature phase. For t — > 0~, the quantities gj approach universal constants that rep- 
resent the low-temperature zero-momentum j-point renormalized coupling constants. A 



simpler parametrization can be obtained if we introduce |86 

M 
M 



M 



so that 



where 



772 

6J r =—^B{u), (5.13) 



w 2 



T^T C - H^O M 2 l; d 

The scaling function B(u) has the following expansion 



w 2 = lim lim ^ . (5-14) 



where 



b(u) = h u -i) 2 + Y: -Mu - iy, (5.i5) 



Vj = JL_. (5.16) 
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B. The four-point zero-momentum renormalized coupling 



The four-point coupling g = g 4 plays an important role in the field-theoretic perturbative 
expansion at fixed dimension [57], which provides an accurate description of the critical 



region in the symmetric phase. In this approach, any universal quantity is obtained from 
a series in powers of g (^-expansion), which is then resummed and evaluated at the fixed- 
point value of g, g* (see e.g. Refs. [^7j.[50] ) . The theory is renormalized at zero momentum 
by requiring 

T^(p) = Z- 1 [M 2 + p 2 + 0(p 4 )] , (5.17) 
r (4) (0,0,0,0) = Z~ 2 Mg. (5.18) 

When M — > the coupling g is driven toward an infrared-stable zero g* of the corresponding 
Callan-Symanzik /^-function 



DM 



(5.19) 

90, A 



In this context a rescaled coupling is usually introduced (see e.g. Ref. 

9 = ^-g. (5.20) 

An important issue in this field-theoretical approach concerns the analytic properties of 
(3(g), that are relevant for the procedure of resummation of the g-expansion. General 
renormalizat ion-group arguments predict a non-analytic behavior of /3(g) at g = g* |]8"T 
One expects a behavior of the form Qj8j| 



(3(g) = -u(g* - g ) + b 1 (g*-g) 2 + ... + c x (g* - g) 1+ ^ + . . . + d x {g* - g)it + . . . (5.21) 

(A = uu and A 2 are scaling correction exponents). In the framework of the 1/N expansion 
of O(N) 4 models, the analysis of the next-to- leading order of the Callan-Symanzik 
/5-function, calculated in Ref. |89| , shows explicitly the presence of confluent singularities of 
the form (|53T| ). 

In the fixed-dimension field-theoretical approach, a precise determination of g* is crucial, 
since the critical exponents are obtained by evaluating appropriate (resummed) anomalous 
dimensions at g*. The resummation of the (7-expansion is usually performed following the Le 
Guillou-Zinn- Justin (LZ) procedure ||7j, which assumes the analyticity of the /3-function. 



The presence of confluent singularities may then cause a slow convergence to the correct 
fixed-point value, leading to an underestimate of the uncertainty derived from stability 
criteria. 

We have computed g* = g\ from our IHT series by calculating the critical limit of the 
quantity g 4 defined in Eq. (|3.12|) . A description of our analysis can be found in App. |A4 . 



The results are reported in Table |V|. We find good agreement among the results of the 
three improved Hamiltonians, that lead to our final estimate: 

g* = 23.49(4), g* = 1.402(2). (5.22) 
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TABLE VI. Results for g\, r 6 , r 8 , r 10 , c 2 and C3 derived from the analysis of the IHT series 
(see App. [5J). The error is reported as a sum of two terms: the first one is related to the spread of 
the approximants; the second one is related to the uncertainty of the value of A4. 

gl r 6 r 8 r w 10 4 c 2 10 4 c 3 

A 6 = 23.499(16+20) 2.051(7+2) 2.23(5+4) -14(4) -3.582(7+6) 0.085(6) 

A 6 = 1 23.491(21+40) 2.050(5+4) 2.23(5+6) -13(5) -3.574(7+20) 0.086(4) 

spin-1 23.487(18+20) 2.046(2+3) 2.34(5+3) -8(25) -3.568(11+4) 0.090(4) 



TABLE VII. Estimates of g* = 3g*/(16ir). (sc) and (bcc) in the HT estimates of Ref. [15] 
denote simple cubic and bcc lattice respectively. For values marked with an asterisk, the error is 
not quoted explicitly in the reference. 



IHT 



HT 



e-exp. 



(isL 1.397M4H 

nj&j 1.391* yr 



d=3 g-exp. 

Br 



MC 



1.39(3) P 
1.408(12; 
1.49(3) 
1.462(12) 



d-exp. 

1.412(14) y : 



ERG 



n 



1.402(2) 



1.408(7) (sc) |1 
1.407(6) Qx;c) 
1.406(9) 
1.414(6) 
1.459(9) 
1.42(9) 




§4. 



i 



1.23(2 
1.72* k 



Table VII presents a selection of estimates of g* obtained by different approaches. The 
HT estimates of Refs. p4]JT5|j36[1 were obtained by using the RT or appropriate biased ap- 
proximants in order to handle the leading confluent correction. The larger result of Ref. 
|9rj| | could be explained by an effect of the scaling corrections. Field-theoretical estimates 



are reasonably consistent, especially those obtained from a constrained analysis of the 0(e 4 ) 
e-expansion JyJ]. In the d=3 g-expansion approach g* is determined from the zero of fl(g) 
after resumming its available 0(g 7 ) series. The results obtained using the LZ resummation 
method JyJ show a slight discrepancy from our IHT estimates. This difference can explain 
the apparent discrepancy found in the determination of 7. Indeed, the sensitivity of 7 to 
g*, quantified in Ref. j|4| through d^/dg* ~ 0.18, tells us that changing the value of g* 
from 1.411 (which is the value obtained from the zero of (3(g)) to 1.402 shifts 7 from 1.2396 
to 1.2380, which is much closer to the IHT estimate 7 = 1.2371(4). Similarly for u, using 
dv/dg* ~ 0.11 |[44|| , v would change from 0.6304 to 0.6294, which is quite acceptable, since a 
residual uncertainty due to the resummation of u(g) is still present. The more general anal- 
ysis of the ^-expansion of Ref. [IE] leads to a smaller value g* = 1.40, with an uncertainty 
estimated by the authors to be about 1%. In Table |V11| we also report estimates obtained 
by approximately solving the exact renormalization group equation ]49j48[] (ERG), and from 
a dimensional expansion of the Green's functions around d = |91| (d-exp.). Concerning 
Monte Carlo (MC) results, we mention that the result of Ref. |90| has been obtained by 



studying the probability distribution of the average magnetization (see also Ref. for a 
work employing a similar approach). The other estimates have been obtained from fits to 
data in the neighborhood of (3 C . In Ref. [IB] Monte Carlo simulations were performed using 
the Hamiltonian (|1.2j ) with A6 = and A4 = 1, which is close to its optimal value. A fit to 
the data of #4, kindly made available to us by the authors, gives the estimate g* = 23.41(24) 
(i.e. g* = 1.397(14)), which is in agreement with our IHT estimate. In Ref. [0 a finite-size 



scaling technique is used to obtain data for large correlation lengths, then the estimate of 
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g\ is extracted by a fit taking into account the leading scaling correction. The Monte Carlo 
estimates of Refs. [94,95] were larger because the effects of scaling corrections were neglected, 

IT41. A more complete list of references regarding this issue can 



as already observed in Ref. 
be found in Ref. |14|]. 



C. Higher-order zero- momentum renormalized couplings 



To compute the HT series of the effective-potential parameters r 2 j defined in Eq. Q5.8|) , 
we rewrite them in terms of the zero-momentum connected 2j-point Green's functions X2j 

as 



rio 



10- 



X6X2 

xi 



280 - 56^^ + XsX2 



xl 



xl 



2 2 2 

15400 - 4620^^ + 126^^ + 120- XsX2 



Xl 



xi 



xl 



XioXl 

xi 



(5.23) 
(5.24) 
(5.25) 



etc... Details of the analysis of the series are reported in App. |A 4j . Combining the results 
reported in Table [VI], we obtain the following estimates: 



r 6 = 2.048(5), 
r 8 = 2.28(8), 
r 10 = -13(4). 

From the results for r 2 j we can obtain estimates of the couplings 



9e 

98 

9w = 



= 9 2 r 6 
= 9 3 r s 
9 4 r w 



1130(5), 
2.96(11) x 10 4 , 
-4.0(1.2) x 10 6 . 



(5.26) 
(5.27) 
(5.28) 



(5.29) 
(5.30) 
(5.31) 



In the literature several approaches have been used for the determination of the couplings 
g 2 j. Table |VIII| presents a review of the available estimates of r&, r§ and rio- We also 
mention the estimate r 10 = —10(2) we will obtain in Sec. |V11| by studying the equation of 
state. The agreement with the field-theoretic calculations based on the e-expansion [44,99 
and on the d 
Ref. 



=3 (7-expansion pi| is good. Precise estimates of r 2 j have also been obtained in 
(see also Ref. ^9|) by ERG, although the estimate of gl by the same method is not 
as good. Additional results have been obtained from HT expansions P^|J9^|J9T[] and Monte 
Carlo simulations [p0| , [9~5]| of the Ising model. The Monte Carlo results do not agree with 
the results of other approaches, especially those of Refs. p2| , |95l , which are obtained using 
finite-size scaling techniques. But one should consider the difficulty of such calculations due 
to the subtractions that must be performed to compute the irreducible correlation functions. 
A more complete list of references regarding this issue can be found in Refs. |23|,p4],|99[ . 
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TABLE VIII. Estimates of Tij- When the original reference reports only estimates of g2j (see 
Refs. [9C,95,9?J), the errors we quote for r2j have been calculated by considering the estimates of 
52 j as uncor related. For values marked with an asterisk, the error is not quoted explicitly in the 
reference. 
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VI. THE TWO-POINT FUNCTION 

The critical behavior of the two-point correlation function G(x) of the order parameter is 
relevant to the description of critical scattering phenomena, which can be observed in many 
experiments, such as light and X-ray scattering in fluids, magnets... In Born approximation 
the cross section for particles of incoming momentum p^ and outgoing momentum pf is 
proportional to the component k = pf — p^ of the Fourier transform of G(x): 

r fi oc G(pf -p^. (6.1) 

As a consequence of the critical behavior of the two-point function G(x) at T c , 

G(k) - j^, (6-2) 

the cross section for k — > (forward scattering) diverges as T — > T c . When strictly at 
criticality, Eq. ( |6.2| ) holds for all fc< A, where A is a generic cut-off related to the microscopic 
structure of the statistical system, e.g. the inverse lattice spacing in the case of lattice models. 
In the vicinity of the critical point, where the relevant correlation length £ is large but finite, 
the behavior ( |6.2|) occurs for A ^> k ^> l/£. At low momentum, k -C l/£, experiments show 
that G(x) is well approximated by a Gaussian (Ornstein-Zernike) behavior, 

+ (63) 

where M ~ l/£ is a mass scale denned at zero momentum (for a general discussion see e.g. 
Ref. ||1Q1|| ). Corrections to Eq. ( |OD are present, and reflect, once more, the non-Gaussian 



nature of the Wilson-Fisher fixed point. The above-mentioned experimental observations, 
confirmed by theoretical studies p9| , |102|| , show that they are small. In the following we will 



improve the determination of the critical two-point function at low-momentum using IHT 
series. 
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In order to study the low-momentum universal critical behavior of the two-point function 
G(x) = (4>(x) (f)(0)) , we consider the scaling function 



g(y)= X /G(k), 



y = k 2 /M 2 



(M = l/£ and £ is the second- moment correlation length) in the critical limit k, M 
with y fixed. The scaling function g(y) can be expanded in powers of y around y = 0: 

oo 

g{y) = l + y + ^CiV*. 



(6.4) 
-> 

(6.5) 



Other important quantities which characterize the low- momentum behavior of g(y) are the 
critical limit of the ratios 



S M = MlJM 1 
S z ee X M 2 /Z t 



gap; 



(6.6) 
(6.7) 



where M gap (the mass gap of the theory) and Z gap determine the long-distance behavior of 
the two-point function: 



G(x) 



H^P_ e - A/gap M. 



47r la: I 



(6.8) 



The critical limits of Sm and Sz are related to the negative zero y of g(y) closest to the 
origin by 



Sm 
Sz 



-yo, 
dg{y) 



dy 



(6.9) 
(6.10) 



y=yo 



The coefficients q can be related to the critical limit of appropriate dimensionless ratios 
of spherical moments of G(x) (as shown explicitly in App. |A 1| ) and can be calculated by 
analyzing the corresponding HT series. Some details of the analysis of our HT series are 
reported in App. [A 4j . In Table |V]] we report the results for A6 = 0, 1 and the spin-1 model. 
We obtain the estimates 



and the bound 



c 2 = -3.576(13) x 1(T 4 , 
c 3 = 0.87(4) x 10" 5 , 

- 10~ 6 < c 4 < 0. 



(6.11) 
(6.12) 



(6.13) 



The constants q and Sm can also be calculated by field-theoretic methods. They have been 
computed to 0(e 3 ) in the framework of the e-expansion [103 ], and to 0(g A ) in the framework 
of the d=3 g-expansion [39|. In Table [IX| we report the results of constrained analyses of the 
0(e 3 ) e-expansion of q and Sm — 1, using exact results in d = 2, 1 (Sm = 1 and c, = in 
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TABLE IX. q and Sm — 1 obtained from 0(e 3 ) series: unconstrained analysis (unc) and 
analyses constrained in dimensions d = 1,2. 



10 



{Sm - 
10 4 c 2 
10 5 c 3 
10 6 c 4 



unc 



d= 1 



-4.4(1.0) 
-4.3(9) 
1.13(27) 
-0.50(13) 



-3.3(8) 
-3.2(8) 
0.84(22) 
-0.37(10) 



d 



-3.3(5) 
-3.3(4) 
0.76(17) 
-0.32(8) 



d= 1,2 



-3.24(36) 
-3.30(21) 
0.69(10) 
-0.27(5) 



TABLE X. Estimates of Sm and Cj. (sc) and (bcc) denote the simple cubic and the 
body-centered cubic lattice respectively. 



C-2 



IHT 



-3.576(13) x 10 



-4 



HT 

-3.0(2) xlO" 4 



e-cxp. 



O.U\^ j A. ±U ||Q^|| 

5.5(1.5) xlO" 4 (sc) M 
7.1(1.5)xlQ- 4 (bcc) |02| 
1.0(1) xlO -5 IE1 



-3.3(2) xlO" 



<i = 3 g-exp. 



-4.0(5) xlO" 



0.87(4) x 10 



-5 



1.0(l)xl0~ 5 
0.5(2) xlO" 5 (sc) flOl 
0.9(3) xl0~ 5 (bcc) fllO 



0.7(1) xl0~ 5 1.3(3) xlO" 



C4 



-10~ 6 < c 4 < 



-0.3(l)xl0- 



-0.6(2)xl0" 



S M 0.999634(4) 



0.99975(10) p 



0.99968(4) 



0.99959(6) 



d — 1; two-dimensional values will be reported in Table [XI]) and following the method of Ref. 
14 1 . Since the constants q are of order 0(e 2 ), we analyzed the O(e) series for Q/e 2 . Errors 
are indicative since the series are short. In Table [X| we compare the estimates obtained by 
various approaches: they all agree within the quoted errors. 
As already observed in Ref. 



39 1, the coefficients show the pattern 



Ci < Cj_i < ... < c 2 < 1 



for 



i > 3. 



(6.14) 



Therefore, a few terms of the expansion of g(y) in powers of y should be a good approximation 
in a relatively large region around y — 0, larger than |y| < 1. This is in agreement with the 
theoretical expectation that the singularity of g(y) nearest to the origin is the three-particle 



cut |[L04| , tL03|| . If this is the case, the convergence radius r g of the Taylor expansion of g(y) 



is r g = 9Sm- Since, as we shall see, Sm ~ 1, at least asymptotically we should have 



1 

— ( 

9 



-Ci. 



(6.15) 



This behavior can be checked explicitly in the large- limit of the N- vector model ||39|| . In 
two dimensions, the critical two-point function can be written in terms of the solutions of a 
Painleve differential equation [|105|| and it can be verified explicitly that r g = 9Sm- In Table 
X| we report the values of Sm and q for the two-dimensional Ising model. 

Assuming the pattern ( |6.14|) , we may estimate Sm and Sz from C2, C3, and C4. Indeed 



from the equation g(yo) = 0, where yo = —Sm, we obtain 
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TABLE XI. Values of Sm and q for the two-dimensional Ising model in the high- and 
low-temperature phase. 



high temperature |10q| low temperature 



S M = 0.999196337056 S M = 0.399623590999 

c 2 = -0.7936796064 x 10~ 3 c 2 = -0.42989191603 

c 3 = 0.109599108 xl0~ 4 c 3 = 0.5256121845 

c 4 = -0.3127446 x 10~ 6 cj = -0.8154613925 

c 5 = 0.126670 xl0~ 7 C5 = 1.422603449 

c 6 = -0.62997 xl0~ 9 Cg = -2.663354573 



S M = l + c 2 -c 3 + c 4 + 2c 2 2 + ... (6.16) 
5 z = l-2c2 + 3c 3 -4c4-2^ + ... (6.17) 



where the ellipses indicate contributions that are negligible with respect to C4. In Ref. |39| 



the relation (|6.16|) has been confirmed by a direct analysis of the HT series of Sm- From 
Eqs. ( |6.16| ) and fl6.17p we obtain Sm = 0.999634(4) (from which we can derive an estimate 
of the ratio = f+ p /f + = 1.000183(2), cf. Eqs. (gJ5|) and (fB6|)) and S z = 1.000741(7). 

We can also use our results to improve the phenomenological model proposed by Bray 
1 103 1 . If we parametrize the large- 2/ behavior of g(y) as [|107|| 



9^ = ^2 ( x + y( i-a)/(2,) + ^m) ' (6 - 18) 

then, by using our estimates of the critical exponents and the phenomenological function of 
Ref. [ 103|| , we obtain the following values for the coefficients: 



A x w 0.918, A 2 w 2.55, A 3 w -3.45. (6.19) 

Estimating reliable errors on these results is practically impossible, since it is difficult to 
assess the systematic error due to the many uncontrolled simplifications that are used. It is 
however reassuring that they are in reasonable agreement the e-expansion predictions 



Ax w 0.92, A 2 w 1.8, A 3 w -2.7, (6.20) 

and with the results of a recent experimental study |DJ 

Ax = 0.915(21), A 2 = 2.05(80), A 3 = -2.95(80). (6.21) 

Bray's phenomenological expression makes also predictions for the coefficients Cj. The pat- 
tern ( |6.15D is built in the approach. We find c 2 = —4.2 • 10~ 4 and c 3 = 1.0 • 10~ 5 , in good 
agreement with our IHT estimates. Therefore, Bray's expression provides a good descrip- 
tion of g(y) for small and large values of y. However, in the intermediate crossover region, 
as already observed in Ref. the agreement is worse: Bray's interpolation is lower by 
20-50% than the experimental result. 

In the low-temperature phase, for y — > 0, the two-point function also admits a regular 
expansion of the form (|6.5| ). However, the deviation from the Gaussian behavior is much 
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larger. The leading coefficient c 2 is larger than c 2 by about two orders of magnitude ||108 
Moreover, by analyzing the low-temperature series published in Ref. [ 109|| one gets = 
0.938(8) (and correspondingly = fg ap /f~ = 1.032(4)). Thus Sj^ shows a much larger 
deviation from one (the Gaussian value) than the corresponding high-temperature phase 
quantity Sm- The two-dimensional Ising model shows even larger deviations from Eq. ( |6.3p . 
as one can see from the values of and a[ reported in Table |X|. Notice that in the 
low-temperature phase of the two-dimensional Ising model the singularity at k 2 
of G(k) is not a simple pole, but a branch point ||105|| . As a consequent 
analytic for \y\ > S^, and therefore the convergence radius of the expansion around y 



—M 2 

gap 

g{y) is not 




is S M . For discussions of the analytic structure of g(y) in the low-temperature phase of the 
three-dimensional Ising model, see e.g. Refs. | ]108| , |103| , |104|J110[ |. 



VII. THE CRITICAL EQUATION OF STATE 

A. The parametric representation 

The critical equation of state provides relations among the thermodynamical quantities 
in the neighborhood of the critical temperature, in both phases. From this equation one 
can then derive all the universal ratios of amplitudes involving quantities defined at zero- 
momentum (i.e. integrated in the volume), such as specific heat, magnetic susceptibility, 
etc... 

From the analysis of IHT series we have obtained the first few non-trivial terms of the 
small-field expansion of the effective potential in the high-temperature phase. This provides 
corresponding information for the equation of state 

H oc t^ s F(z), (7.1) 

where z oc Mt~^ and, using Eq. ( |5.9| ), 

F(z) = ^ = z + |z 3 + £ F 2m+1 z 2m+1 (7.2) 



with 



2m- 1 



(2m - 1)! 



r 2m - (7.3) 



The function H(M,t) representing the external field in the critical equation of state (|7.1|) 
satisfies Griffith's analyticity: it is regular at M = for t > fixed and at t = for M > 
fixed. The first region corresponds to small z in Eq. ( [7.1[ ), while the second is related to 
large z, where F(z) can be expanded in the form 

F(z) =z s J2 F™z- n/p . (7.4) 

n=0 

Of course F°° are universal constants. 
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To reach the coexistence curve, i.e. t < and H = 0, one should perform an analytic 
continuation in the complex t-plane Jl],|23| . The spontaneous magnetization is related to the 
complex zero Zq of F(z). Therefore the description of the coexistence curve is related to 
the behavior of F(z) in the neighbourhood of z . In order to obtain a representation of the 
critical equation of state that is valid in the whole critical region, one may use parametric 
representations, which implement in a simple way all scaling and analytic properties. One 
parametrizes M and t in terms of R and 9 



M = m R l3 9, 
t = R(l-9 2 ), 

H = h R^ 5 h{9), (7.5) 

where ho and mo are normalization constants. The function h(9) is odd and regular at 9 = 1 
and at 9 = 0. The constant h can be chosen so that h(9) = 9 + 0(9 3 ). The zero of h(9), 
9 > 1, represents the coexistence curve H = 0, T < T c . The parametric representation 
satisfies the requirements of regularity of the equation of state. One expects at most an 
essential singularity on the coexistence curve [ [111 



The relation between h(9) and F(z) is given by 

z = p9(l-9 2 y\ (7.6) 
h(9) = p- 1 (l-9 2 Y S F(z(9)), (7.7) 

9 > 0, and hyperscaling implies that (58 = /3 + j. Notice that this mapping is invertible only 
in the region 9 < 9i, where 9i = (1 — 2 / 3)~ 1 / 2 is the solution of the equation z'{9) = 0. Thus 
the values of 9 that are relevant for the critical equation of state, i.e. < 9 < 9q, must be 
smaller than 9\. This fact will not be a real limitation for us, since the range of values of 9 
involved in our calculations (which will be < 9 2 < 9q < 1.40) will be always far from the 
limiting value 9f ~ 2.88. 

As a consequence of Eqs. (|7.5| ), ( |7.6|) , and ( [7.7|) , we easily obtain the relationships 

(7 ' 8) 

We can therefore treat p as a free parameter, and the scaling relations between physical 
variables will not depend on p, provided that itlq and ho are rescaled with p. In the exact 
parametric equation the value of p may be chosen arbitrarily but, as we shall see, when 
adopting an approximation procedure the dependence on p is not eliminated, and it may 
become important to choose the value of this parameter properly in order to optimize the 
approximation. 

From #0 one can obtain the universal rescaled spontaneous magnetization [p3|l , i.e. the 
complex zero z of F(z), 

z = \z \e- mfi , \zq\ = P 9 (91 (7.9) 

From the function h(9) one can calculate the universal ratios of amplitudes. In App. [B] we 
report the definitions of the universal ratios of amplitudes that have been introduced in the 
literature, and the corresponding expressions in terms of h{9). 
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Expanding h{9) in (odd) powers of 9, 

h(9) = 9 + "£h 2n+1 9 2n+1 , (7.10) 



n=l 



and using Eq. (|7.7Q , one can find the relations among h 2n +i and the coefficients F 2m +i of 
the expansion of F(z). The procedure is explained in App. £], and the general result is: 

n 

h 2 n+l — 2^ Cn^mp m F 2m+ i, (7.11) 
m=0 

where 

i n—m 

J] (2/3m- 7 + k- 1); (7.12) 



(n-m)! 

notice that c nj „ = 1. In general h 2n+ \ depends on 7, /3, and on the coefficients F 2m+ i with 
m < n. 

We shall need the explicit form of the first two coefficients: 

/i 3 = |p 2 -7, (7.13) 
h 5 = §7(7 - 1) + |(2/3 - 7 )p 2 + F 5 p 4 . (7.14) 



B. Approximation scheme based on stationarity 

In Ref. [^| Guida and Zinn- Justin use the first few coefficients of the small- 2 expansion 



of F(z) to get polynomial approximations of h(9) which should provide a description that is 
reliable in the whole critical region. The approximations considered are truncations of the 
small-# expansion of h(p,9), i.e. 



t-i 



h^(p,9) = 9 + Y / h 2n+ i(p)9 2n+1 , (7.15) 



71=1 



where h 2n+ i(p) are given by Eq. ( [7.1 We follow a similar strategy, with a significant 
difference in the procedure adopted in order to fix the value of p . 

By Eqs. (|7.11|) and ( [7.1 2D , the coefficients h 2n+ x(p) included in Eq. (|7.15|) are written in 



terms of the t parameters 7, (3, F5 ... F 2t -\. In practice only the first coefficients of the 
small-6 1 expansion of h(9) are well determined, since we have good estimates only for the 
first few F 2m+ i. Once the order of the truncation has been decided, one may exploit the 
freedom of choosing p to optimize the approximation of h(9). In this way one may hope 
to obtain a good approximation even for small values of t. Ref. [^3| proposes to determine 
the optimal value of p by minimizing the absolute value of h 2t _i(p), i.e. the coefficient of 
the highest-order term considered. The idea underlying this procedure is to increase the 
importance of small powers of 9. Our approach is different. 

Our starting point is the independence on p of the scaling function F(z) and, as a 
consequence, of all universal ratios of amplitudes that can be extracted from it. Of course, 
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this property does not hold anymore when we start from a truncated function (p, 
if we compute universal quantities from a function (p, z) defined by 



FV{p,z)=FV{p,9{p,z)), 



where 



F«(p, 



9 2 ) 



its 



i.e. 



(7.16) 



(7.17) 



and 6(p, z) is obtained by inverting Eq. ( |7.6|) . 

In order to optimize p for a given truncation h®(p,6), we propose a procedure based 
on the physical requirement of minimal dependence on p of the resulting universal function 



F®(p,z). This can be obtained by assuming p to depend on z, i.e. 
requiring the functional stationarity condition 



p = p 



(*) ( 



5FW(p®,z) 
5pW 







and by 



(7.18) 



(see Ref. [H5| and references therein for a similar technique applied to the resummation of 
perturbative power expansions). The non-trivial fact, even surprising at first sight, is that 
the solution p^'(z) of Eq. ( |7.18| ) is constant. In other words, for any t there exists a solution 
pt independent of z that satisfies the global stationarity condition 



dF®{p,. 



dp 



(7.19) 



p=pt 



This is equivalent to the fact that, for any universal ratio of amplitudes R, its approximation 
Ry'(p) (obtained from Fw(p,z)) satisfies the stationarity condition 



dR®{p) 



dp 



0. 



^7.20) 



p=pt 



The proof of Eq. ( |7.19| ) is given in App. y, where we show that the global stationarity 
condition amounts to requiring p t to be a solution of the algebraic equation 



(2(3 - l)p^ - 2 7 + 2t - 2 



(7.21) 



The idea behind our scheme of approximation is that, for any truncation, the stationarity 
condition enforces the physical request that the universal ratios of amplitudes be minimally 
dependent on p. To check the convergence of the approximation, one can repeat the compu- 
tation of universal ratios of amplitudes from the truncated function h^'(pt,9) for different 
values of t, as long as one has a reliable estimate of i^t-i- We have no a priori argument 
in favor of a fast convergence in t of the universal ratios of amplitudes derived by this 
procedure towards their exact values. However, we may appreciate that its lowest-order 
implementation, corresponding to t=2 in Eq. ( |7.15| ), reproduces the well-known formulae of 
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Refs. p4[-[2~6|j, which give an effective optimization of the linear parametric model. Indeed 



we obtain from Eqs. ( [7.21 ) and ( 7.131) the t = 2 solution 



In this case the critical equation of state and all critical amplitudes turn out to be expressible 
simply in terms of the critical exponents /3 and 7. In particular, we found a closed-form 
expression for all -F 2 m+i coefficients (see App. y for a derivation): 

= ^T^^ UW™ — j — k). (7.23) 



mi P2 



k=l 



Wallace and Zia |27| already noticed that the minimum condition of Refs. P^-|2q| was 
equivalent to a condition of global stationarity for the linear parametric model. We have 
shown that such a global stationarity can be extended to other parametric models, regardless 
the linearity constraint, and can be used to improve the approximation. 

The next truncation, corresponding to t = 3, can also be treated analytically. Since it 
sensibly improves the linear parametric model in the 3d Ising case, we shall present here a 
few details. By applying the stationarity condition (|7.21 ) to Eq. (|7.14 ), we obtain 



( T -2/3)(l-7 + 2/3) 



1 



\ 12(4/3-7)^5 I \ (7-2/3)2(1-7 + 2/3) 



x _ 72(2- 7)7(7 -1)(4^- 7)F 5 _ ( ^ 4) 



Universal ratios of amplitudes may be evaluated in terms of they will now depend only on 
the parameters /3, 7 and F§. Notice that the predictions of the t = 2 and t = 3 models differ 
from each other only proportionally to the difference between the "experimental" value of 
F 5 and the value predicted according to Eq. ( [7.231) 

5 " 72 7(7-1)- (T 5) 

(2) 

If we replace F 5 with F 5 in the t = 3 model results, all the linear parametric model results 
are automatically reproduced. In the 3d Ising model, the two values differ by ~ 6%, and 
thus we expect comparable discrepancies for all universal ratios of amplitudes. This can be 
verified from the numerical results that we will present in Sec. [VII D| (see Table [XI I|) . All 
universal ratios of amplitudes obtained from the t = 2 truncation (i.e. the linear parametric 
model), using our estimates of 7 and /3, differ at most by a few per cent from previously 
available estimates. The t = 3 and higher-order approximations are consistent with the 
latter. The apparent convergence in t of the results provides a further important support to 
this scheme. 

It is worth noticing that the parametric representation of the equation of state induces 
parametric forms for such thermodynamics functions as the free energy and the susceptibil- 
ity, as discussed in detail in App. [B 2j . When we assume a truncated form of the parametric 
equation of state, in general only the corresponding free energy function will admit a poly- 
nomial representation. A peculiar and possibly unique feature of our scheme is the induced 
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truncation of the function related to the susceptibility, which turns out to be an even poly- 
nomial of degree 2t in the variable 9. App. |C| contains a more extended discussion of these 
and other properties of the approximation scheme based on the stationarity condition. 

We have introduced our parametric representation assuming independent knowledge of 
F 5 , F 2 t-i- It should be noticed that our results for t = 3 can also be used as a phe- 
nomenological parametrization, fitting the value of F§ on any known universal quantity. As 
we will show in Sec. |V11 Dj , the difference with the linear parametric model of Refs. [33-21) 
is not negligible. On the other hand, our numerical estimates for t — 4 show that the differ- 
ence from t = 3 is too small (compared with both theoretical and experimental precision) 
to justify the introduction of an additional phenomenological parameter F 7 . 



C. e-expansion of the parametric representation 

It is interesting to compare our results with the analysis of the parametric equation of 
state which can be performed in the context of the e-expansion, generalizing results presented 
in Refs. WW^i- 



According to Ref. [p7[ , within the e-expansion it is possible to choose a value po such 
that, for all n > 2, 



^2n+l(Po) =0(e 



(7.26) 



The calculation shows that p — v2- We proved in App. that Eq. ( [7.261 ) keeps holding for 
all choices of p that satisfy the relation p = p + 0(e). We can now e-expand our globally 
stationary solutions for arbitrary t, obtaining 



lim p t 



(7.27) 



As a consequence, any truncation satisfying the stationarity condition is an accurate de- 
scription of the e-expanded parametric equation of state up to O(e') included. 

As a byproduct, we may extract from the linear model relation ( |7.23|) , expanded to 0(e 2 ), 
the coefficients of the e-expansion for i*2m+i) for rn > 2: 



2m+l 



^ ] fmk£ ■ 



(7.28) 



fc=i 



We easily obtained from Eq. (|7.23|) the closed form results 

f-l) m 1 



fml 
fm2 



6m(m — 1) 2 

fml 



1\ m ~ ^ 1 



17 m / m ± * ^ x 
27 ~ ~2 ~ V 3~ + 6/ ^ k 



(7.29) 
(7.30) 



reproducing known results ||23[1 . More generally, knowing the expansion of the coefficients 
-p2m+i to 0(e t ) for m < t is enough to reconstruct all F 2m +i for m > t to the same accuracy. 
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TABLE XII. Universal ratios of amplitudes obtained by taking different approximations of 
the parametric function h(8). Numbers marked with an asterisk are inputs, not predictions. The 



values jO m t are obtained as in Ref. [23|, see text for details. 





h {2> (p2,d) 


7 / r,\ 

h^'(p 3 ,9) 




^ (5, (P5,f) 




h (4) (Pm,4,0) 


p 


1.7358(12) 


1.7407(14) 


1.7289(83) 


1.686(51) 


1.6889(26) 


1.651(30) 


«§ 


1.3606(11) 


1.3879(29) 


1.372(12) 


1.325(53) 


1.3310(13) 


1.295(27) 


Fg° 


0.03280(14) 


0.03382(18) 


0.03374(21) 


0.03366(26) 


0.03378(18) 


0.03370(23) 


\M 


2.825(12) 


2.7937(17) 


2.7970(33) 


2.8012(72) 


2.7955(15) 


2.7992(45) 




0.5222(16) 


0.5316(21) 


0.5295(29) 


0.5261(60) 


0.5303(19) 


0.5276(39) 


u 2 


4.826(11) 


4.752(15) 


4.769(22) 


4.797(47) 


4.764(13) 


4.786(30) 




-9.737(41) 


-8.918(83) 


-9.10(18) 


-9.42(48) 


-9.061(67) 


-9.31(28) 


RJ 


0.05538(13) 


0.05681(16) 


0.05644(32) 


0.0558(10) 


0.05656(14) 


0.05606(53) 


Rc 


0.021976(16) 


0.022488(30) 


0.02235(11) 


0.02211(36) 


0.022387(18) 


0.02220(19) 




7.9789(64) 


7.804(10) 


7.823(18) 


7.847(40) 


7.8146(85) 


7.836(25) 




92.10(23) 


93.91(20) 


93.25(45) 


91.9(1.7) 


93.25(21) 


92.27(83) 


1>3 


6.0116(79) 


6.0561(68) 


6.041(11) 


6.010(39) 


6.0412(73) 


6.018(20) 


V4, 


16.320(55) 


16.121(66) 


16.21(11) 


16.41(24) 


16.239(55) 


16.38(15) 




1.6775(19) 


1.6588(27) 


1.6624(44) 


1.668(10) 


1.6611(25) 


1.6656(60) 


U 2 Rt 


38.505(82) 


37.09(15) 


37.31(26) 


37.64(55) 


37.23(13) 


37.50(35) 




0.4419(13) 


0.4434(13) 


0.4416(18) 


0.4377(56) 


0.4420(13) 


0.4392(29) 




1.9389(48) 


*2.048(5) 


*2.048(5) 


*2.048(5) 


*2.048(5) 


*2.048(5) 


rs 


2.507(31) 


2.402(39) 


*2.28(8) 


*2.28(8) 


2.365(43) 


*2.28(8) 


no 


-12.612(41) 


-12.146(60) 


-10.0(1.5) 


*-13(4) 


-11.80(10) 


-10.98(86) 



D. Results 



As input parameters for the determination of the functions (p, 9) we use the results 
of the IHT expansion: 7 = 1.2371(4), v = 0.63002(23), r 6 = 2.048(5), r 8 = 2.28(8), 
r 10 = -13(4). 

In Table [XII| we report the universal ratios of amplitude as derived from truncations 
corresponding to t — 2,3,4,5. We use the standard notation for the ratios of amplitudes 
(see e.g. Ref. [|113j| ); all definitions can be found in Table XIII . For comparison, we also 
report, for t = 3,4, the results obtained using the procedure of Ref. P3L fixing p to the 
value p m) t which minimizes the absolute value of the O(0 2t ~ l ) coefficient^. Such results are 
very close to those derived from the stationarity condition; this is easily explained by the 
fact that the values of p mjt are close to p t . The errors reported in Table [Xll| are related to the 
uncertainty of the corresponding input parameters (considering them as independent). The 
results for t = 2,3,4 suggest a good convergence and give a good support to our analysis. 
The results for t = 5, although perfectly consistent, are less useful to check convergence, due 
to the large uncertainty of Fg. In Table |X111| we report our final estimates, obtained using 
h^(p4,6); all the approximations reported in Table [XI I| are consistent with them, except 
t = 2. 

We should say that the method of Guida and Zinn- Justin to determine the optimal p 
leads to equivalent results, and shows an apparent good convergence as well. However, we 
believe that the global stationarity represents a more physical requirement, and it is more 



4 As already noted in Ref. [23|, for t 
never reaches zero. 



4 the minimum of hi{p) is zero, while, for t = 3, h§(p) 
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TABLE XIII. Summary of the results obtained in this paper by our high-temperature calcula- 
tions (IHT), by using the parametric representation of the equation of state (IHT-PR), by analyzing 
the low-temperature expansion (LT), and by combining the two approaches (IHT-PR+LT). Nota- 
tions are explained in App. |H. 



IHT 



IHT-PR 



LT 



IHT-PR+LT 



7 
v 

a 

i) 

P 

6 

a 

>'6 

rs 

rio 

U = 

U 2 = 

u 4 = 

R~ = 

i?4 = 
R 3 = 
i? 4 = 



A ' A 

c + /c- 
ct /c A 

aA+C+/B 2 
aA~C-/B 2 
-C+B 2 /{C + f = 
v 3 = -C 3 ~B/(C- 
C A B 2 /(C~f 

— i?4 + 3i?| 



\ z o\ 



Qi 5 = R X = C+B s - l /{5C c ) & 



5 4 = 
w 2 = 

% = 
Q + ~= 
Rt=- 

Q~ = 
Q c = 

93 = 
9 A = 

Qt = 
Qs = 

Q2 = 



1/3 



cf. Eq. q[A 

g = -C 4 + /[(C+) 2 (/+) 3 ] 
C-/[B 2 (f-f] 

f + /f- = (™ 2 U 2 Rt/g+ 
= aA+(f+) s = R+R+/gt 
e (Q+) 1/3 

e aA-(f-) 3 = R-/w 2 
B 2 (f+)yC+ = Q+/R+ = Ri/gj 
wv 3 
w 2 va 

:f+ lf+ 
- J gap / J 

~ f gap / f 

= f + If 
— J gap/ J i 

/gap/ f 

i (/7/ + ) 2 ""c+/c- 



gap/ J gap 
c c 



1.2371(4) 
0.63002(23) 
0.1099(7) 
0.0364(4) 
0.32648(18) 
4.7893(22) 
0.0208(12) 
2.048(5) 
2.28(8) 
-13(4) 



23.49(4) 



-10(2) 
0.530(3) 
4.77(2) 

-9.1(2) 
0.0564(3) 
0.02235(11) 
7.82(2) 
6.041(11) 

93.3(5) 
16.21(11) 
1.662(5) 
0.0337(2) 



4.75(4) 



0.01880(8) 
0.2659(4) 

0.3330(10) 



1.961(7) 



0.00471(5) 



1.000183(2) 



1.032(4) p 



13.17(6) 
77.0(8) 



1.901(10) 

1.024(4) 

1.195(10) 
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amenable to a theoretical analysis of its convergence properties. Moreover, as we have shown, 
it has the linear parametric model of Refs. p4|-p7H as the lowest-order approximation. 

Estimates of other universal ratios of amplitudes can be obtained by supplementing 
the above results with the estimates of w 2 = C~ /[B 2 (f~) 3 } and = f gap /f~ obtained 
by an analysis of the corresponding low-temperature expansion. The results so obtained 

The low-temperature expansion of w 2 can be 

The 



are denoted by IHT-PR+LT in Table [Xffl 
calculated to 0{u 21 



on the cubic lattice using the series published in Refs. ||109|,|114 



results reported in Table |XII1| were obtained by using the Roskies transform in order to 
reduce the systematic effects due to confluent singularities jjTJJ]. 

We also consider a parametric representation of the correlation length. Following Ref. 
we write 
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CJx = R- 



"gap 



(9). 



We consider the simplest polynomial approximation to a{9) and a gap (9): 

a(6) « a (l + c6 2 ) , 

a gap($) ~ a gap,0 (l + Cgap^" 



(7.31) 
(7.32) 



(7.33) 
(7.34) 



where the constants c and c gap can be determined by fitting the quantities U$ and . 
Then, using Eqs. ( |7.31| ), ( |7.32| ), and the parametric representation of the equation of state, 
one can estimate the universal ratios of amplitudes and Q2 defined in Table XIII . Notice 



that, given the equation of state, the normalization ao is not arbitrary, but it may be fixed 
using the zero-momentum four-point coupling g±: 



a = (h /p) 1/3 (m /p)- 5/3 (gl)- 2 / 3 



-5/3 



(7.35) 

where h , m and p have been introduced in Eqs. ( |7.5| ) and (|7.6|). Notice that a depends 
only on the ratios h /p and m /p, as it is required of a physical quantity. Moreover one has 
Ogap.o = (Q^) 2 a o- In order to check the results obtained from the approximate expressions 
( [T.33| ) and (|7.34|) , we also considered the following parametric representation [ 102 1 : 



-2 



R 2v b(9), 



(7.36) 
(7.37) 



and the corresponding polynomial approximations truncated to second order. The results 
for and Q2 obtained by this second representation are perfectly consistent with those 
from the first one. Our final estimates of and Q2 derived by the above method are 
reported in Table [Xlll . 

In Table |XIV| we compare our results with other approaches. We find a good overall 
agreement. Our results appear to substantially improve the estimates of most of the univer- 
sal ratios considered. In Table PCIV| we have collected results obtained by high-temperature 
and low-temperature expansions (HT,LT), Monte Carlo simulations (MC), field-theoretical 
methods such as e-expansion and various kinds of expansions at fixed dimension d = 3, 
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TABLE XIV. Estimates of the quantities in Table Kill] by various approaches. The experi 



mental data are taken from Ref. [55], unless otherwise stated, ms denotes a magnetic system; bm a 
binary mixture; lv a liquid-vapor transition in a simple fluid. For values marked with an asterisk, 
the error is not quoted explicitly in the reference. 



IHT-PR 



HT,LT 



MC 



e-exp. 



d=3 exp. 



experiments 



U 



0.530(3) 



0.523(9 
0.51 



0.560(10) 
0.550(12) 
0.567(16) 



= 

10 



0.527(37) 
0.524(10) 



3^ 

118 119 



0.537(19) 
0.540(11) 
0.541(14) 



[L2X1 

[I2i| 



bm 
lv 
ms 
lv (56 



4.95(15 r Upi 
5.01* $0 : 4 



(pip 



lv 



0.56(2) 
0.50(3) 
0.51(3) 
0.53 



+ 8 



57 0.538(17) 



U 2 



4.77(2) 



4.75(3) [rMT 
4.72(11)1 123 1 



4.73 
4.9* 
4.8* 



llfi 



fj|||l25 



4.79(10) 
4.77(30) 
4.72(17) 



[L2E| 



bm 
lv 
ms 
ms 



4.4(4) 
4.9(2) 
5.1(6) 



u 4 


-9.1(2) 


-9.0(3) \3t\^ 


-8.6(1.5) 


11 


-9.1(6) \i4 


\ 




Rt 


0.0564(3) 


0.0581(10 


) | 


6j 0.0569(35) JlJ 


0.0574 
0.0594 


,20) [ 


B 

121 


bm 0.050(15) 
lv 0.047(10) 


*t 


7.82(2) 


7.94(12) 


m 


8.24(34) 




7.84* | 




Ri = v- A 


6.041(11) 


6.44(30) 


113 


96[ 5.99(5) | 






6.08(6) 




il 







R~ 93.3(5) 107(13) 


9d 


113 






v 4 16.21(11) 15.8(1.4) 8 


6l 


Qi" 1.662(5) 1.57(23) 


127 113| 1.648f36l 

1.67* [US 


M 1.669 18) Q bm 1.75C301 
LIS | 1.7* [L21I lv 1.69(14) 



4.71(5) 




1.961(7) 



1.96(1U1 
1.96* [102 



1.95(2) 
2.06(1) 



22 

[129 



0.0193(10) [@ 



1.91* [[121] 
0.0197* Jl3lj|ll9 



2.013(28) [ 126 ] bm 



1.93(7) 
1.92(15) 



0.01880(8) 



0.0202(9) p 
0.01880(15) [116 



0.01968(15) |13 



56 



0.0174(32) 
0.023(4) 



0.00471(5) 0.00477(20 



113 



0.0463(17)J 10 



0.3330(10) 0.324(6) 



0.328(5) 



0.331(9) [h21| 



bin 
lv 



0.33(5) 
0.35(4) 



g~ 13.17(6) 


13.9(4) pe| 


13.6(5 


M 


13.06( 


12) |86| 


94 77.0(8) 


85* [36| 


108(7) 


133| 


75(7) [ 


86| 



1.000183(2) 1.0001* 



1.00016(2) 



1.00021(3) 



1.024(4) 



1.007(3) [11 



.031(6) |ll0||l3l| 



1.032(4) 
1.037(3), 



1.195(10) 



1.17(2) [1127 pl3 



1.13* 



and experiments. Concerning the HT,LT estimates, we mention the recent Ref. [|113j| where 
a review of such results is presented. The agreement with the most recent Monte Carlo 
simulations is good, especially with the results reported in Ref. 
cise. However we note that the estimates of Uq reported in Ref. 



122 



which are quite pre- 
are slightly larger. 

Moreover there is an apparent discrepancy with the estimate of of Ref. [ |133|| . It is worth 
mentioning that the result of Ref. [|123jl was obtained simulating a four- dimensional SU(2) 
lattice gauge model at finite temperature, whose phase transition is expected to be in the 
3d Ising universality class. Field-theoretical estimates are in general less precise, although 
perfectly consistent. We mention that the results denoted by ll d = 3 exp." are obtained from 
different kinds of expansions: g-expansion [^5] , ^ , p!2l|J132|| , minimal renormalization without 
e-expansion ||120| , |135|| , expansion in the coupling u = 3w 2 defined in the low-temperature 
phase [126]. In Refs. [^,0 Guida and Zinn- Justin used the d = 3 g- and e-expansion to 



calculate the small-field expansion of the effective potential and a parametric representation 
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of the critical equation of state. We also mention the results (not included in this table) of an 
approach based on approximate solution of exact renormalization equations (see e.g. Refs. 
p|,p|p6|). Some results can be found in Ref. (T3§: U 2 S 4.29, g% S 15.24, Q{ s = 1.61 
and = 1.86. 

We report in Table [XiV| experimental results for three interesting physical systems ex- 
hibiting a critical point belonging to the 3d Ising universality class: binary mixtures, liquid- 
vapor transitions and uni-axial antiferromagnetic systems. A review of experimental data 
can be found in Ref. [55]. Most of the results shown in Table |XIV| were reported in Refs. 



P^| , |122| | . They should give an overview of the level of precision reached by experiments. 

For sake of comparison, in Table [XV| we report the universal ratios of amplitudes for the 
two-dimensional Ising model. The purely thermal results are taken from Ref. ||105|| , where 
the exact two-point function has been written in terms of the solution of a Painleve equation. 
Q + and Q~ have been computed by us solving numerically the differential equations reported 
in Ref. [ 105|| . The ratios involving amplitudes along the critical isotherm can be obtained 
using the results reported in Ref. ||138|| . For the quantities that are not known exactly, we 
report estimates derived from the high- and low-temperature expansion. Such estimates are 
quite accurate and should be reliable because the leading correction to scaling is analytic, 
since the subleading exponent A is expected to be larger than one (see e.g. Ref. ||139[ | and 

[U! 



references therein). In particular the available exact calculations [ |105|| for the square-lattice 
Ising model near criticality have shown only analytic corrections to the leading power law. 
Therefore the traditional methods of series analysis should work well. 
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APPENDIX A: GENERATION AND ANALYSIS OF THE 
HIGH-TEMPERATURE EXPANSION FOR IMPROVED HAMILTONIANS 



1. Definitions 



Before discussing the series computation, let us define all the quantities we are interested 
in and fix the notation. 

Starting from the two-point function G(x) = ((f)(0) <p(x)), we define its spherical moments 



m 2j = £(x 2 ) J G(x) 



(Al) 



(x = mo) an d the first non-spherical moments 



X 



2 V 



z 4 -§(* 2 ) 2 



G(x) 



(A2) 
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TABLE XV. Universal ratios of amplitudes for the two-dimensional Ising model. Since the spe- 
cific heat diverges logarithmically in the two-dimensional Ising model, the specific heat amplitudes 
A^ are defined by Ch ~ —A ^ Int. 

7 7/4 

v 1 

U = A+/A- 1 

U 2 = C+/C- 37.69365201 

R+ = A+C+/B 2 0.31856939 

R~ = A-C-/B 2 0.00845154 

Q{ s = R X = C+B s - l /{5C c ) s 6.77828502 

w 2 = C-/[B 2 (f-) 2 } 0.53152607 

U{ = f + /f~ 3.16249504 

^€gap = /gap/ /gap ^ 

Q+ = 4+(/+) 2 0.15902704 

Q- = A-{f-) 2 0.015900517 

Qt = fg^/f + 1.000402074 

Q\ = /|a P // c 1.0786828 

Q% = f gap /r 1.581883299 

Q2 = {f c /f + ) 2 - f) C+/C c 2.8355305 



9t = 9 = -C+/[(C+) 2 (/+) 2 ] 14.694(2) @0 

r 6 3.678(2) 

r 8 26.0(2) H 

Ho 275(15) 



v 3 = -C^B/(C-) 2 33.011(6) m 

v 4 = -C4 B 2 /(C-) 3 + 3u| 48.6(1.2 



(where x™ = X^scf). 

The second-moment correlation length is denned by 

£2 = M -2 = (A3) 

The coefficients c, of the low- momentum expansion of the function g(y) introduced in Sec. |VT| 
can be related to the critical limit of appropriate dimensionless ratios of spherical moments, 
or of the corresponding weighted moments m2j = iri2j/x- Introducing the quantities 

^ = ^jTT)! M2iM2i ' (A4) 

one can define combinations of u 2 j (that we will still call q to avoid introducing new symbols) 
having; critical limit: 

c 2 = l-w 4 + ^M 2 , (A5) 
c 3 = 1 - 2u A + u 6 - ^M 4 , (A6) 
c 4 = 1 - 3w 4 + u\ + 2u e -u 8 + ^M 4 + okoM 6 , (A7) 
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etc... Notice that the terms proportional to powers of M 2 do not contribute in the critical 
limit M — > 0, but they allow us to define improved estimators [|3^| . Indeed in the lattice 
Gaussian limit, defined by the two-point function 

G(k) = ¥ ^- 2 , P = £4sin 2 (^/2), (Ag) 

Q = independently of M, and not only in the critical limit M — > 0. 
The zero-momentum connected Green's functions are defined by 

X2j= E (<P{0)<p{x 2 )...<p{x2j-i)<P{x2 3 ))c] (A9) 

X2,—,X2j 

in particular, X2 = X- 



2. Linked cluster expansion 

We computed the high-temperature expansion by the linked cluster expansion (LCE) 



technique. A general introduction to the LCE can be found in Ref. ||140|| . We modeled the 



application of the LCE to 0(AQ-symmetric models after Ref. [141 . 

In order to perform the LCE for the most general model described by the Hamiltonian 
( |1.2| ), we parametrize the potential V(<j) 2 ) in terms of the "single-site moments" | 141 | 



m 2k = ,) 2 1 N N ~ 1+2k , J k = I dxx k exp(-V(x 2 )). (AlO) 



2 k T^N + k) Jn-i 

We compute our series for fixed N, leaving all wi2k as free parameters: each term of the 
series is a polynomial in m2k with rational coefficients. 

With the aim of computing as many terms of the series as possible, we adopted all the 
technical developments of Ref. [ I4"2 | , and we introduced more improvements of our own; in 



this Appendix, we will only describe these; readers not familiar with technical details of the 
LCE should consult Refs. [|141|J142j| . 



As discussed in Ref. ||142|| , Sec. 3, the LCE requires a unique representation of graphs; it 
is convenient to implement this by defining a canonical form for the incidence matrix. The 
reduction to canonical form of a graph with V vertices requires in principle the comparison of 
the V\ incidence matrices obtained by permutation of vertices, which is clearly unmanageable 



for large graphs; even with the introduction of the "extended vertex ordering" of Ref. ||142 
this operation remains the dominant factor in the computation time; therefore, we devoted 
a large effort to the optimization of this aspect of the computation. On one hand, we have 
perfected the extended vertex ordering, and we are able to recognize inequivalent vertices 
much more often. On the other hand, we search for (a subgroup of) the symmetry group of 
the incidence matrix, which allows us to perform explicitly only one permutation for each 
equivalence class. Altogether, the largest sets of vertices which are explicitly permuted are 
of size 5 or less (except for a few hundred diagrams requiring 6, and a handful requiring 
permutations of 7 or 8 elements). 
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The next most computer-intensive operation is the computation of embedding numbers 
and color factors; it is optimized by "remembering" each computed value in a table, compat- 
ibly with available memory. This is crucial for color factors, which are computed recursively, 
and very effective for embedding numbers. 

The problem of handling integer and rational quantities which do not fit into machine 
precision is solved by using the GNU multiprecision (gmp) library. Neither multiprecision 
nor polynomials in m2k are necessary for the most expensive sections of the computation; 
therefore they have little impact on the computation time. 

In order to speed up the handling of search and insertion into ordered sets of data, we 
make extensive use of AVL trees (height-balanced binary trees) (cfr. e.g. Ref. ||143|| , Chapt. 
6.2.3). AVL trees are used to manipulate graph sets, multivariate polynomials, and tables 
of embedding numbers and of color factors. 

The LCE is dramatically simplified by restricting actual computations to the set of 
one-particle irreducible graphs. One must however establish the relationship between the 
usual moments and susceptibilities and their irreducible parts. To this purpose, we found it 
convenient to define a generating functional of irreducible moments (irreducible momentum- 
space 2-point functions): 

G lpi (p) = ]Texp(zp ■ x) G lpi (x), (All) 

X 

where G lpi (x) is the irreducible-graph contribution to the field-field correlation. One may 
then prove the relationship 

[G(p)]- 1 = [G^ipT 1 - 2 / 3]Tcosp J . (A12) 

i=l 

By expanding both sides of Eq. ( |A12|) in powers of pf, it is trivial to establish all desired 
relationships, both for spherical and for non-spherical moments. 

We have calculated \ and m 2 to 20th order, and the other moments of the two-point 
functions to 19th order. We have calculated Xa to 18th order, xe to 17th order, xs to 16th 
order, and xio to 15th order. Using Eqs. ( |3.12|) , ( |5.23|) , (|5.24|) , and (|5.25|) , one can obtain 



the HT series corresponding to zero-momentum four-point coupling g± and the quantities 
r2j we have introduced to parametrize the effective potential. 
It is useful to factorize out the leading dependence on (3: 

n 

= /7$>#; (A13) 

i=0 

the values of r and n are summarized in Table pCVJ| . In the following we will analyze the 
series normalized to start with 0(/3°), i.e. ao + /3a± + ... 

We have checked our series for X2n and 777,2 against the available series of the standard 



Ising model (see e.g. Refs. [T^j3g]); in this special case our only new result is the 18th order 



coefficient of the expansion of Xi- 



171450770247965944104542584 
ai8(X4j - 32564156625 ' ^ ' 
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TABLE XVI. Summary of normalization and length of our IHT series. 



O r n 

~X 20 

1 19 

9a -| 17 

r 6 17 

r 8 16 

Ho 15 

c 2 4 13 

c 3 3 13 

c 4 2 13 



We have checked the (new) series for m 4 and m 6 by changing variables to v — tanh (3 and 
verifying that all coefficients become integer numbers. 

It would be pointless to present here the full results for an arbitrary potential: the 
resulting expressions are only fit for further computer manipulation. For the three potentials 
we are interested in, we computed m>2k by numerical integration (to 32 digit precision or 
higher). The coefficients a« of the HT series for A6 = and A4 = 1.1 and for A6 = 1 and 
A4 = 1.9 are reported in Tables [XV11| and |X VlTT] respectively. The series for the spin-1 model 
defined by Eq. (|1.4j ), with D = 0.641, are reported in Table [XIX . 



3. Critical exponents 

In order to determine the critical exponent 7 from the nth order series of % (n — 20 in 
our case) we used quasi-diagonal first, second and third order integral approximants (IAl's, 
IA2's and IA3's respectively). 

IAl's are solutions of the first-order linear differential equation 

Pi(x)f(x) + P (x)f(x) + R(x) = 0, (A15) 

where the functions Pi(x) and R(x) are polynomials that are determined by the known nth 
order small- a; expansion of f{x). We considered [mi/mo/fc] IAl's with 

mi + m + k + 2 > n — p, 

Max[[(n -p- 2)/3j - q, 2] < mi,m Q ,k< \(n-p-2)/3] +q, (A16) 

where m 1; m , k are the orders of the polynomial Pi, Pq and R respectively. The parameter 
q determines the degree of off-diagonality allowed. Since the best approximants are expected 
to be those diagonal or quasi-diagonal, we considered sets of approximants corresponding to 
q = 3. For a given integer number p, only approximants using n terms with n > n > n — p 
are selected by (|A16| ). In our analysis we considered the values p — 0, 1. 



IA2's are solutions of the second-order linear differential equation 

P2(x)f"(x) + Px(x)f(x) + P (x)f(x) + R(x) = 0. (A17) 
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TABLE XVII. IHT series for A 4 = 1.1, A 6 = 0, in the notation of Eq. (|A13|) . 



1 


X 


m2 


JTJ4 


me 





0.5308447611308816674 











1 


1.6907769625206168952 


1.6907769625206168952 


1.6907769625206168952 


1.6907769625206168952 


2 


4.8619910172171602715 


10.770481113538254496 


28.721282969435345322 


86. 163848908306035967 


3 


13.927143445449825611 


48.231864289335148202 


219.75546850876176116 


1168.8527451895890195 


■1 


38.903779467013842036 


188.05200645833592066 


1229 . 89 1 2709035843796 


9662.7116928978662527 


5 


108.53608309082300208 


675.87104272683156690 


5832.6111092485591617 


61434.000941321774177 


6 


299.26769419241406575 


2309.0702869335329515 


24922.163253442953617 


332091.05589205489865 


7 


824.59140866646260666 


7605.2927810004027352 


99136.392920539417813 


1607441.8553355166225 


8 


2256.9464691160956608 


24394. 253637363046531 


374199.26312351789406 


7181423.6410152125984 


9 


6174.4168460205032479 


76627.209216394784414 


1356978.6661855495156 


30189095.008610024020 


10 


16819.879385593690953 


236799.07667260657765 


4767392.0099162457155 


120979394.99942695905 


11 


45803 222727034040360 


721928 56294418965891 


16324368 332600723787 


466443860 83571669282 


12 


124363.26835432977776 


2176629.7412550662641 


54721422.909210908059 


1742031434.2158738860 


13 


337573.74963124787949 


6500509.9999124014627 


180180889.59217314660 


6334414835.2189134577 


14 


914347.51881247417342 


19258195.825678901432 


584286774.86828839115 


22515209536.395375901 


15 


2476042.0363235919004 


56653368.820414865245 


1869882856.0802601457 


78474598792.492084546 


16 


6694111.3933182426254 


165647959.63348214176 


5915646880.4853995649 


268883672555.38373998 


17 


18094604.163418613844 


481713425.26839816320 


18526280962.623843359 


907568914837.87389317 


18 


48847832.893538572297 


1394159442.3995129568 


57500007423.420861038 


3022858145406.1028395 


19 


131848611.02423050678 


4017559436.5586218326 


177034534120. 13444060 


9949496764882.0179674 


20 


355511932.47075480765 


1 1532862706. 754267638 






i 


X4 


X6 


X8 


Xio 





—0.3285640660980093563 


1.0162413020868264428 


—6.8394260743547676540 


79.348906365011205720 


1 


—4.1859963164157358115 


25.898050097125655128 


—286.46285381596139097 


4926.0498824266722100 


2 


—30.741724996686361292 


334.99400966656728040 


—5647.7911946128472205 


136420. 74020456605398 


3 


— 176.53992137927974557 


3085.9534765513322166 


—75010.539622248383529 


2448228.6617818973385 


1 


—873.19795113342604018 


22962.676420970454675 


—772240.19869798230118 


33066114.972427672856 


5 


-3914.7539681033628945 


147391.51074600633293 


—6637987.3112014615524 


364115447.75387624436 


6 


— 16340.897140343714854 


848040.97611583569625 


—49819302.751247368336 


3432806859.5766444770 


7 


—64653.470284596876205 


4484314.0227663949556 


—336218289.35412445864 


28624021619.321695067 


8 


—245234.96688614286659 


22168058.403461618364 


—2082668486.5197392127 


215985650524. 12674877 


9 


-899257.84241359785957 


103719954.32133484536 


-12019618575.926765056 


1499822923734.9247399 


10 


-3206654. 1029660245146 


463530204.01159412237 


-65362330202.998798196 


9707913378518.1866160 


11 


-11170819.137408319309 


1992634695.4221763266 


-337846725947.59204461 


59157394345479.522579 


12 


-38148679.051940544866 


8285182135.0991054145 


-1671339685352.6240713 


342085313031 1 14.78303 


13 


- 128071730. 18983471414 


33466763808.875455279 


-7957523515361.0083551 


1889271366312617.3703 


11 


-423602975.57444761466 


131799942528.27775884 


-36629909644439.962983 


10018120044594804.748 


15 


- 1382909952.8269485885 


507558776082.59290672 


-163635571877509.11631 


51230076380872446.008 


16 


-4462746050.5347000940 


1915992506452.6475673 


-711670345253605.66031 




17 


-14253923929.146690502 


7104558940304.9779228 






18 


-45107295178.923296542 









We considered [r^/mi/mo/A:] IA2's with 
m 2 + m i + m o + k + A > n — p, 

Max[[(n-p-4)/4j -q,2]< m 2j m 1} m ,k < [(n-p-4)/4] + q, (A18) 

where m 2 ,mi,m , k are the orders of the polynomial P%, P\, Pq and R respectively. Again 
the parameter q determines the degree of off-diagonality allowed, and we consider sets of 
approximants corresponding to q — 2. 

IA3's are solutions of the third-order linear differential equation 

P 3 (x)f"'(x) + P 2 {x)f"{x) + P 1 {x)f{x) + P (x)f(x) + R(x) = 0. (A19) 

We considered [7713/7^2/ wii/mo/ A;] IA3's with 

m 3 + m 2 + rrii + m o + k + 6 > n — p, 

Max [\_(n — p — 6)/5j — q, 2] < m 3 , m 2 , mi, mo, k < \(n — p — 6)/5] + q. (A20) 
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TABLE XVIII. IHT series for A 4 = 1.9, A 6 = 1, in the notation of Eq. (|A13|) 



l 


X 


m2 


JTJ4 


me 





0.4655662671465330507 











1 


1.3005116946285418792 


1.3005116946285418792 


1.3005116946285418792 


1.3005116946285418792 


2 


3.2867727127185864813 


7. 2656925005834656853 


19.375180001555908494 


58. 125540004667725483 


3 


8.2759367806706747301 


28.571904795057895903 


130.05174486699400177 


691.57352659837378755 


■1 


20.312529070243777317 


97.872800052605172528 


638.95885613696143199 


5016. 1872942503349323 


5 


49.794483203882242272 


309.06347996559814197 


2661.2142140348739306 


27994. 182028926323346 


6 


120.62506581503040408 


927.75225884375746610 


9988.4118850888212117 


132873.57808445335775 


7 


292.01512515747113432 


2684.8185172507997345 


34904.461272349558702 


564842.72150668529468 


8 


702. 16648194884216675 


7566.3506306935978232 


115746.99820735924226 


2216511.1678444057633 


9 


1687.6428772292658916 


20882.094002198352809 


368765.50096064887178 


8184889.2748500168622 


10 


4038.7968304000570954 


56696.707453278695709 


1138240.4288922759401 


28813837.691500628877 


11 


9662 2729936573649464 


1 51 863 269821 82684635 


■H424267 0781090006844 


9759591 3830641 471 31 


12 


23047. 113364799224307 


402271.59171959125975 


10084793.611992545198 


320214529.50558023009 


13 


54959.286082391555970 


1055489.7357595885681 


29174001.507844433518 


1022942183.3309332937 


14 


130774.53710516575464 


2747201.8540323841754 


83116951.940376136421 


3194364390.8578671172 


15 


311110.31236133784486 


7100121.4902629495116 


233696696.59215911285 


9781470144. 1262582176 


16 


738903.79800753751923 


18238478.826399953589 


649550962.77873855718 


29444730149.930149105 


17 


1754633.3818777485475 


46596282.414619468624 


1787187905.0885303101 


87315557903.862730209 


18 


4161219.5431591307916 


118476705.93261076006 


4873253901.0331542961 


255504682223.76959562 


19 


9867152.2571694621736 


299943150.21579862108 


13181855628.291951390 


738841191899.23449896 


20 


23372660.203375142541 


756429264.37368967452 






i 


X4 


X6 


X8 


XlO 





-0.2477796481363361878 


0.6474598784982264035 


-3.7217906767025756640 


37.013405312193356664 


1 


—2.7685883005851708908 


14.535362477178409902 


— 137.07512921231932765 


2014.8224919025960131 


2 


— 17.877560582981223820 


165.70820207079800267 


—2380.8220701405889901 


49093.219551382446508 


3 


—90.322805654272940178 


1345.3663489182523069 


—27871.670743361006737 


776123.99513244411557 


1 


—393.02065426992610848 


8820.7362245901247130 


—252935.23058597306885 


9238516.9347113237042 


5 


- 1549.8790048885508754 


49871.943864090224636 


-1916250.8876237147160 


89671059.994338257912 


6 


—5689.6695707018226942 


252682.48564411688205 


— 12673128.761295226766 


745143260.15631620580 


7 


— 19795.290447377859541 


1176295.0381597851008 


—75349625.332312049631 


5475799689.5566930641 


8 


—66017.086587684140210 


5118124.4639866386541 


—411110139.91029709454 


36408447077.502609221 


9 


-212822.79826473961810 


21072833.518843848267 


-2089380048.3140434305 


222743859156.05535943 


10 


-667124.90202267593913 


82859426.257839273325 


-10003697691.834647146 


1270011540667.3854746 


11 


-2042817.2956666846385 


313350700.53757124062 


-45518090240.849158501 


6816107886501.6215858 


12 


-6131769.6806946007770 


1146009747.0203113140 


-198194900838.32752342 


34708877742973 .050540 


13 


- 18092493.782248051554 


4071313753.6304055865 


-830439274834.30354664 


168778204360271.32357 


11 


-52592383.442563299219 


14100289090.570869350 


-3363665155462.8191575 


787889916506321.20775 


15 


-150889630.61165148373 


47747821664.976757633 


-13220617132779.560155 


3546559365689180.8887 


16 


-427911822.28715837917 


158482857336.16489875 


-50583038860319.049950 




17 


-1201047651.7022285867 


516675010346.30966717 






18 


-3339910306.5273359851 









We considered sets of approximants with q = 2. 

As estimate of (3 C and 7 from each set of IA's, we took the average of the values corre- 
sponding to all non-defective approximants listed above. The error we quote is the standard 
deviation. Approximants are considered defective when they present spurious singularities 
close to the real axis for Re (3 % (3 C . More precisely we considered defective the approximants 
with spurious singularities in the rectangle 

x min < Rez < x max , \lmz\ < y mauX (A21) 

where z = (3/(3 c . The special values of x min , x max and y max are fixed essentially by stability 
criteria, and may differ in the various analysis. One should always check that, within a 
reasonable and rather wide range of values, the results depend very little on the values 
of x min , x max , and y max . The condition ( |A21| ) cannot be too strict, otherwise only few 
approximants are left. In this case the analysis would be less robust and therefore less 
reliable. In the analysis of the critical exponents we fixed x m ; n = 0.5, a; max = 1.5 and 
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TABLE XIX. IHT series for the spin-1 model at D = 0.641, in the notation of Eq. ( |A13|) . 



1 


X 


mi 


7714 


7776 





0.5130338416658140921 











1 


1.5792223361663016211 


1.5792223361663016211 


1.5792223361663016211 


1.5792223361663016211 


2 


4.4354864269385182067 


9.7223340236143130564 


25.926224062971501484 


77.778672188914504451 


3 


12.419491590434972647 


42.346809834993589030 


191.98340105778667094 


1019.9725391572417242 


■1 


33.791261055416279831 


160.98579902461367310 


1043.4445460560890431 


8166.3561104899782974 


5 


91.841711514904201711 


564.00666776024395394 


4815.3070263215122748 


50394.717412597175181 


6 


246.45252550898370464 


1877.9601316327756011 


20036.716263907461491 


264816. 649443 14665340 


7 


661.02891687705666766 


6026.3998915743002481 


77637.764412960183453 


1247156.1972913768310 


8 


1760.2485042097784968 


18829.052334178552758 


285474.20654850401533 


5423869.9597014869631 


9 


4685.9023162918237050 


57600.716340522523798 


1008432.9506590528288 


22201293.939743630952 


10 


12417.403239568002677 


173324.21633191405061 


3450900.8811683331178 


86641778.653884117489 


11 


32897 780017066145091 


^144^0 09422679607285 


1 1 50871 2 41 8722055563 


3253321 57 496850481 74 


12 


86884.987650751268743 


1509919.6022647962218 


37570398.804037029475 


1183310452. 1936542300 


13 


229424.66023796502223 


4389267. 7826683406013 


120463565.05045829114 


4190401181.7453112798 


14 


604434. 19630362500628 


12656102.176997207352 


380359458.94391796446 


14504931859.074136829 


15 


1592166.9939650416411 


36234119.703394609364 


1185138499.8785854689 


49231114687.307644210 


16 


4186778.2089201191352 


103100174. 19850927716 


3650157580.5925683557 


164257098550.90866386 


17 


11008100.036706697979 


291755921.33251277078 


11128154630.539032580 


539841179613.24402666 


18 


28904025.427069749972 


821639224.91959656471 


33620312241.342947618 


1750686524137.2717781 


19 


75884596.302083003892 


2303832207.2589619187 


100754983048.94595013 


5610143692354.6457755 


20 


199011100.35574405792 


6434727599.1288912159 






i 


X4 


X6 


X8 


XlO 





—0.2765773264173367185 


0.6159505110893571460 


—2.9991913236448423261 


25.639540414004216349 


1 


—3.4054446789491065718 


15.965825500576724169 


— 131.09722874868945636 


1673.3057480542363748 


2 


—24.570159964990933901 


209. 17712832988124972 


—2673.1181642081994199 


48487.693538927944393 


3 


— 139.04496175368819963 


1941.9314061871553985 


—36418.744278635469309 


903571.05878070985095 


1 


—676.69722203050291907 


14481.161465713991310 


—381670.70958440882592 


12572456.840428388939 


5 


-2980.2850555040833295 


92716.594093707793992 


-3318236.6956045245656 


141625856.39478899275 


6 


— 12201. 170328998898251 


530071.30088368454939 


—25054434.002489226202 


1357653370. 1441467390 


7 


—47288.798001661947371 


2776591.5196504503986 


— 169364792.32951428214 


11451941832.665412740 


8 


— 175521.23292975181789 


13563137.626324266820 


— 1047051748.2193837515 


87035901992.179110029 


9 


-629296.76861756838811 


62579486.626554750976 


-6012877960.8045829709 


606512902870.52199480 


10 


-2192523.1579334804197 


275333754.56012219662 


-32454249850.587780052 


3927249840355.5223369 


11 


-7458569.3308517637566 


1163638939.7290930525 


-166147643973.94469398 


23876098793430 . 398990 


12 


-24861185.382016374834 


4751112394.9532757501 


-812614050679.29384404 


137426012654323.52351 


13 


-81432521.735317351753 


18826990715.719661290 


-3819167751639.5573175 


753928052790611.90171 


11 


-262699929.49158922336 


72675832066.845002904 


-17330670355620.376973 


3964188356224407.7875 


15 


-836235938.80657455801 


274128146369.26786810 


-76232334051550.698537 


20070120847424617.843 


16 


-2630659142.3082426455 


1012932475761 . 7342902 


-326121067148253.92416 




17 


-8189047799.9075180970 


3674558299661.5473376 






18 


-25252383565.446882139 









Umax = 0.5. Sometimes we also eliminated seemingly good approximants whose results were 
very far from the average of the other approximants (more than three standard deviations). 

As a further check of our analysis we used the fact that \ must present an antiferro- 
magnetic singularity at [3f = —(3 C with exponent 1 — a [j37| , cf. Eq. (|4~1~1) . We verified the 



existence of a singularity at (3 ~ — f3 c and calculated the associated exponent. In some 
analyses we selected the approximants with a pair of singularities (3f and (3 C such that 
Pc + Pf — £ Pc, and extracted the estimates of (3 C) 7 and 7 af from them. As in Ref. 0, we 
also considered IA's where the polynomial associated with the highest derivative of f(x) is 
even, i.e. it is a polynomial in x 2 . We will denote them by b a fIA's. This ensures that if x c 
is a singularity of an approximant then — x c will also be a singular point. 

In Table [XX| we present the results for some values of the parameters p, q, e introduced 
above (when the value of e is not explicitly shown it means that the corresponding constraint 
was not implemented). We quote the "ratio of approximants" r-app (I — s)/t, where t is the 
total number of approximants in the given set, I is the number of non-defective approximants 



43 



(passing the test ( |A21| )), and s is the number of seemingly good approximants which are 
excluded because their results are very far from the other approximants; / — s is the number 
of "good" approximants used in the analysis; notice that s <C I, and I — s is never too small. 
We found the IA2 analysis to give the most stable results, especially with respect to the 
change of the number of terms of the series considered. Therefore we consider the IA2 results 
to be the most reliable. Moreover, IAl's are not completely satisfactory in reproducing the 
antiferromagnetic singularity when its presence is not biased. In the biased analyses where 
C is forced, IAl's, IA2's and IA3's give almost equivalent results. 

The results are quite stable, and the value of 7 af is always consistent with 1 — a ~ 0.89. 
From the results of Table [XXL combining the results of the IA2 and IA3 analyses (selecting 



the results denoted by IA2 g=2i p=o, b a fIA2g =2iP =o, IA3g=2, p =o and b a fIA3 g= 2, p =o) we obtain the 
following estimates 

P c = 0.3750973(14), 7 = 1.23732(24) for A 6 = 0, A 4 = 1.10, (A22) 
P c = 0.4269780(18), 7 = 1.23712(26) for A 6 = 1, A 4 = 1.90, (A23) 

and 

P c = 0.3856688(20), 7 = 1.23680(30) for spin-1, D = 0.641. (A24) 
Taking also into account the uncertainty of A4 and D*, we arrive at the estimates of Table 



III. Notice that the value of p c at A4 = 1.10 and A6 = is in agreement with the Monte 



Carlo estimate of Ref. f20|| , i.e. /3 C = 0.3750966(4) (where according to the author the error 



does not include possible systematic errors). From the analysis of the antiferromagnetic 
singularity using the b a fIA's we obtain the following estimate for a: 

a = 0.105(10), (A25) 



which is in good agreement with the much more precise estimate ( |4.5|) obtained assuming 
hyperscaling. 

In order to determine v from the analysis of the HT series of £ 2 , we followed the suggestion 
of Ref. PS|, i.e. to use the estimate of /3 C derived from the analysis of x i n order to bias 
the analysis of the series of £ 2 . We analyzed the 19th order series of £ 2 / 'P and employed 
biased integral approximants (bg c IA). For instance, biased IA2's can be obtained from the 
solutions of the equation 

(1 - x/p e ) P 2 (x)f"(x) + P 1 (x)f'(x) + P (x)f(x) + R(x) = 0. (A26) 

In this case we considered the approximants satisfying the conditions 

m 2 + mi + mo + k + 3 > n — P, 

Max[[(n -p- 3)/4j - q, 2] < m 2 ,m 1 ,m ,k < [(n-p-3)/4] + q, (A27) 

where, as before, rrii and k are the order of the polynomials Pi and R respectively. We also 
tried doubly-biased IA2 (b± ( a c IA2) where also a singularity at — p c is forced using solutions 
of the equation 
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TABLE XX. Results of various analyses of the 20th order IHT series for x- r-app is explained 
in the text. In the biased analyses forcing the value of (3 C the error is reported as a sum: the first 
term is related to the spread of the approximants at the central value of f3 c , while the second one 
is related to the uncertainty of the value of f3 c and it is estimated by varying (3 C . 



approx 



r-app 



/3c 



7 



7 



A 6 = 0, A 4 = 1.08 



A 6 = 0, A 4 = 1.10 



IA2 9 =2,p=o 
IA3,j=2,p=o 
b a f 1^3^=2,^=0 

IAl g =3, p =0 

IA1 



IA1 



g=3,p=0,e=10- 2 
=10- 3 




q—3,p— 0,e 
b afIAl 9 =3 iP = 

IA2 5 =2,p=0 

IA2 



'g=2,p=0,£=10- 
<j=2,p=0,e=10- 



A 6 = 0, A 4 = 1.12 



IA2, 

b afIA2 9 =2,p= 

IA2g = 2,p=l 

IA2 o=2,p=l,e 
IA3 g =2,p=0 
IA3 q=2,p=0,e 

b af IA3 9 =2 lP = 

IA3 q = 2,p=l 

b )3 McIAl g= 3 
b )3 McIA2 9 =2 
b ± «McIA2 9= 
bj_^McIA3g- 
IA2 g =2,p=o 
IA3 g =2,p=o 

t>afIA3 g =2.p= 





= 10- 3 

io- 3 



p=0 
p=0 
2,p=0 
2,p=0 



(78 - 2)/85 
(62 - l)/65 
(63 - l)/65 

30/37 
(26 - l)/37 

1/37 
(24 - 2)/37 
(78 - 2)/85 
(69 - l)/85 

16/85 
(54 - l)/85 
(151 - 3)/165 

25/165 
(62 - l)/65 

21/65 
(63 - l)/65 
(99 - 3)/100 

38/48 
(114 - 3)/115 
(90 - 2)/100 
(61 - l)/63 
(77 - 2)/85 
(62 - l)/65 
(62 - l)/65 



0.3760701 
0.3760699 
0.3760703 
0.3750945 
0.3750955 
0.3750956 
0.3750937 
0.3750975 
0.3750974 
0.3750975 
0.3750989 
0.3750974 
0.3750975 
0.3750971 
0.3750983 
0.3750976 
0.3750957 
0.3750966 
0.3750966 
0.3750966 
0.3750966 
0.3741203 
0.3741199 
0.3741203 



22) 
21) 
18) 
27) 
18) 

46) 
18) 
21) 
19) 
31) 
26) 
16) 
21) 
12) 
18) 
55) 

4) ii 

4) [2C] 
4) [2C] 
4) |C] 
16) 
21) 
17) 



1.23717(36) 

1.23713(37) 

1.23719(30) 

1.23684(42) 

1.23699(26) 

1.23701 

1.23673(69) 

1.23734(29) 

1.23733(35) 

1.23735(31) 

1.23753(46) 

1.23731(39) 

1.23735(25) 

1.23728(36) 

1.23749(19) 

1.23734(30) 

1.23704(98) 

1.23718(7+7) 

1.23720(2+7) 

1.23719(3+7) 

1.23719(8+7) 

1.23748(26) 

1.23742(36) 

1.23747(29) 



-0.376(6) 
-0.377(4) 

-0.375(3) 

-0.3760(8) 

-0.3754 

-0.376(6) 
-0.375(2) 
-0.3751(2) 

-0.376(4) 
-0.3751(2) 
-0.375(4) 
-0.3750(2) 

-0.375(3) 



-0.374(4) 
-0.375(4) 



-1.0(7) 

-1.0(3) 

-0.903(10) 

-0.9(3) 

-0.77(8) 

-0.81 

-0.887(11) 

-1.0(6) 

-0.9(2) 

-0.90(4) 

-0.906(18) 

-1.0(5) 

-0.90(3) 

-1.0(5) 

-0.92(2) 

-0.902(11) 

-1.0(4) 



-1.0(6) 
-1.0(3) 
-0.904(9) 



A 6 = 1, A 4 = 1.86 
A 6 = 1, A 4 = 1.90 



IA2 9 =2,p=o 

IAl 9 =3,p=o 
bafIAl,j=3,p = 
IA2,j=2,p=0 







IA2 



q=2,p=0,e=10- 3 



A 6 = 1, A 4 = 1.94 



b a f IA2q=2, p = 

IA2g = 2,p=l 

IA2 a=2,p=l,e 
IA3,j=2,p=0 
kafIA3,j=2,p= 
IA2q = 2,p=Q 





= 10- 3 




(82 - 3)/85 

37/37 
(29 - l)/37 
(83 - 2)/85 

23/85 
(69 - 2)/85 
(156 - 4)/165 

31/165 
(63 - l)/65 
(63 - l)/65 
(82 - 3)/85 



0.4307605 
0.4269723 
0.426972(7) 
0.4269779 
0.4269779 
0.4269791 
0.4269777 
0.4269777 
0.4269782 
0.4269786 
0.4232606 



22) 
81) 
) 

24) 
32) 
32) 
32) 
30) 
24) 
25) 
21) 



1.23676(30) 

1.2363(11) 

1.2363(9) 

1.23711(34) 

1.23710(45) 

1.23725(48) 

1.23707(42) 

1.23708(41) 

1.23714(38) 

1.23719(36) 

1.23738(30) 



-0.431(5) 
-0.429(7) 

-0.428(5) 
-0.4271(2) 

-0.428(5) 

-0.4271(2) 

-0.429(6) 

-0.423(6) 



-1.1(5) 

-0.7(5) 

-0.890(16) 

-1.0(5) 

-0.90(3) 

-0.903(16) 

-0.9(5) 

-0.89(4) 

-1.0(4) 

-0.906(10) 

-1-1(9) 



spin-1, D = 0.633 
spin-1, D = 0.641 



spin-1, D = 0.649 



IA2 ? =2,p=o 

bafIA2,j = 2,p = 
bafIAl 9 =3,p= 

IA2,j=2,p=o 
IA2, 
b af IA2 
IA2 g =2,p=i 

IA3q=2,p=0 

b at IA3 9 =2,p= 
IA2,j=2,p=o 
b a f IA2q = 2,p= 



-o=2,p=0,e 
2q=2 



=10~ 3 




(82 - 4)/85 
(73 - 4)/85 
(30 - 3)/37 
(73 - 2)/85 
(37 - 2)/85 
(73 - 4)/85 
(146 - 3)/165 
(62 - 5)/65 
(61 - 2)/65 
(82 - 3)/85 
(73 - 4)/85 



0.3845065 
0.3845076 
0.3856634 
0.3856681 
0.3856685 
0.3856691 
0.3856669 
0.3856686 
0.3856682 
0.3868365 
0.3868377 



27) 
17) 

41) 
33) 
17) 
16) 
53) 
36) 
38) 
32) 
16) 



1.23683(34) 

1.23698(25) 

1.2360(6) 

1.23674(38) 

1.23678(23) 

1.23687(25) 

1.23655(69) 

1.23678(58) 

1.23673(60) 

1.23660(36) 

1.23676(24) 



-0.384(8) 



-0.386(3) 
-0.3860(3) 

-0.385(10) 
-0.385(4) 

-0.386(7) 



-1.1(1.2) 

-0.905(12) 

-0.887(11) 

-0.9(4) 

-0.87(3) 

-0.905(12) 

-1.2(2.0) 

-1.1(5) 

-0.911(13) 

-1.2(1.0) 

-0.905(12) 
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TABLE XXL 


Results of the analysis of the 19th order IHT series for £ 2 /(3. r- 


-app is explained 


in the text. The error 


is reported as a sum of two terms. The first term is related to the spread of 


the approximants 


at the central value of f3 c , 


while the second one 


is related to the uncertainty of 


the value of j3 c and it 


is estimated by varyinj 


? 0c- 










approx 


r-app 




7 


\ — n \ a — 1 ns 

y\Q — U, — l.UO 




b ( 3 c IA2 9=2 ,p=o 


(54 - 2)/70 


0.63004(2+11) 




\ „ _ n X/i — 1 1 n 




b^ c IAl g =3 lP =o 


35/37 


0.63012(2+10) 








b | g c IA2 9= 2 ! p=o 


51/70 


0.63016(3+9) 








t>±/3 c IA2 9= 2,p=0 


(53 - 3)/55 


0.63015(5+10) 


-0.88(5) 






b / 3 c IA3 g= 2,p=o 


28/35 


0.63009(14+10) 








b/3 c IA2 g= 2 l p = i 


(86 - 4)/132 


0.63016(3+9) 








b^McIA2 9 =2,p=o 


(47 - 2)/55 


0.63012(3+3) 




A 6 = 0, A 4 = 1.12 




b /3t ,IA2 9= 2 i p=o 


(49 - 2) /70 


0.63027(3+9) 




A 6 = 1, A 4 = 1.86 




hrj TA9 i a 

u p c lrt -^g— l,p— 


("fin — 1 "I /70 


629781' 1 +1 1 1 




A 6 = 1, A 4 = 1.90 




b/3 c IAl q= 3 iP= o 


34/37 


0.63000(2+12) 








b / 3 c IA2 9= 2,p=o 


(61 - 2)/70 


0.63003(2+11) 








b±/3 c IA2 9= 2,p=o 


(55 - 2)/55 


0.63003(7+11) 


-0.87(16) 






b J 3 c IA2 9= 2,p=i 


(113 - 4)/132 


0.63003(3+10) 








b/3 c IA3 9= 2, p =o 


(27 - l)/34 


0.62988(17+14) 




A 6 = 1, A 4 = 1.94 




b ( 3 c IA2 9= 2,p=o 


(55 - 2)/70 


0.63023(2+11) 




spin - 1, D = 0.633 




b / 3 c IA2 9= 2,p=o 


(67 - l)/70 


0.62998(2+13) 




spin - 1, D = 0.641 




b/3 c IAl 9 =3,p=0 


33/37 


0.62988(5+15) 








b^ c IA2 9= 2,p=o 


(66 - l)/70 


0.62990(2+13) 








b^ c IA3 9= 2,p=o 


24/35 


0.62981(12+19) 




spin -!,£> = 0.649 




b | g c IA2 9= 2 1 p=0 


(66 - l)/70 


0.62982(2+13) 





(l - x 2 /ff) P 2 (x)f"(x) + P 1 {x)f{x) + P (x)f(x) + R(x) = 0. (A28) 

In Table [XXI| we report the results of some of the analyses we performed. In the case of 
the b+ / 3 c IA2 analysis we also report the exponent at the antiferromagnetic singularity which 
turned out to be always consistent with 1 — a. The error of v is given as a sum of two 
terms: the first one is computed from the spread of the approximants at /3 C , the second one 
is related to the uncertainty of [3 C . We quote as our final estimates: 

v = 0.63015(12) for A 6 = 0, A 4 = 1.10, (A29) 



v = 0.63003(13) for A 6 = 1, A 4 = 1.90. (A30) 

and 

v = 0.62990(15) for spin-1, D = 0.641. (A31) 

Taking also into account the uncertainty of A 4 we arrive at the estimates of Table ffTTL We 
mention that unbiased IA analyses of the 19th series of £ 2 //3 give consistent but less precise 
estimates of f3 c (cf. Eqs. (|A22|) and flA2p ) and v. 

As a check of our results, we performed a biased analysis of \ an d £ 2 at A6 = and A 4 = 
1.10, using the value f3 c = 0.3750966(4) obtained in Ref. |2(| by Monte Carlo simulations 



based on finite-size scaling techniques. Although the author of Ref. [20] says that the error 
on f3 c does not include systematic errors, we used it as a check and found (see Tables |XX| 
and |XXT1 ) 7 = 1.23720(2 + 7) and v = 0.63012(3 + 3) (the first error is related to the spread 



of the approximants at f3 c = 0.3750966 and the second one to the error on /? c ), which are 
perfectly consistent with our final estimates reported in Table fTJ. 
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TABLE XXII. Results for rj obtained using the CPRM: (a) applied to £ 2 and \ (20 orders); 



(6) applied to £ 2 //3 and x (19 orders), r-app is explained in the text. 





approx 


r-app 




AQ — U, A4 — l.Uo 


(a~\ KT A 9 

^a; Dli\Zq = 2,p=0 


^yo — ± ) / 1 10 


fl 0997zl7l 
U.UZZ 1 <±^±o J 




KT A 9 rt 
l^Dj DLI\Zq — 2,p— 


00/ 1 u 


n c\ooqa( fi\ 
u.uzzy^o ) 


AQ — U, A4 — ±.1U 


kt a 9 ^, 

l^aj Dl/\Zq — 2,p— 


yo/ 110 


fl 099&IV 1 






^lOu — 1J/ too 


n n99sn/'i fi^ 






1^4/ — 1 J/Dl 


n n99fini'Q'7 s \ 




fK\ KT A 1 

^DJ DlAlq = 3 )P =0 


9S /'37 
Zo/ O 1 


n n9Qnri('s\ 




^K^ KTA9 n 


36/70 


n D9^ni i^k 1 ! 

U.UZOU-L^O J 




(hS KTA9 1 


/I ^9 


n n9miQH 9 s ! 




yuj uirvoq — 2,p— 


31/34 


02*31 1 f22^ 


\ ,. — \ , — 1 1 9 


fnS KTA9 ^ 
l^a; uii-iZg— 2,p— 


97 /115 


fl 0998"^ 1 ^ 




fK\ KTA9 ^ n 


3^ /7fl 


n n9^n^^Q^ 
u.uzouo^y i 




f«\ KT A 9 ,i 
[A) DLJ\Z q — 2,p— 


/Q/1 C\ /I 1 C 

^y^± — oj/ 110 


n n99R7/'i 9 s ! 




(b) bIA2g= 2 ,p=o 


39/70 


0.02285(12) 


A 6 = 1, A 4 = 1.90 


(a) bIA2g^2,p-o 


(90-2)/115 


0.02278(12) 




(b) bIA2g^ 2 ,p-o 


37/70 


0.02298(11) 


A 6 = 1, A 4 = 1-94 


(a) bIA2g= 2 ,p=o 


(92-2)/115 


0.02288(13) 




(b) bIA2g= 2 ,p=o 


32/70 


0.02312(10) 


spin-1, D = 0.633 


(a) bIA2g= 2 ,p=o 


(84- 1)/115 


0.02292(40) 




(b) bIA2g^ 2 ,p-o 


36/70 


0.02316(22) 


spin-1, D = 0.641 


(a) bIA2g^ 2 ,p-o 


(85-2)/115 


0.02288(40) 




(b) bIA2g= 2 ,p=o 


37/70 


0.02312(22) 


spin-1, D = 0.649 


(a) bIA2 g =2,p=o 


(84- 1)/115 


0.02285(43) 




(b) bIA2g= 2 ,p=o 


37/70 


0.02307(22) 



In order to obtain an estimate of 77 without using the scaling relation 7 = (2 — r))v, 
we employed the so-called critical point renormalization method (CPRM). The idea of the 
CPRM is that from two series D(x) and E(x) which are singular at the same point Xq, 
D(x) = J^idiX 1 ~ (xq — x)~ 5 and E(x) = Yl,i^i^ % ~ (^0 — x )~ e , one constructs a new series 
F(x) = Ysi{di/ Gi)x % ■ The function F(x) is singular at x = 1 and for x — > 1 behaves as 
F(x) ~ (1 — x)~^, where <fi = 1 + 8 — e. Therefore the analysis of F(x) provides an unbiased 
estimate of the difference between the critical exponents of the two functions D(x) and E(x). 
The series F(x) may be analyzed by employing biased approximants with a singularity at 
x c — 1. In order to check for possible systematic errors, we applied the CPRM to the series 
of £ 2 / '(3 and \ (analyzing the corresponding 19th order series) and to the series of £ 2 and x 
(analyzing the corresponding 20th order series). We used IA's biased at x c — 1. In Table 
pCXITj we present the results of the analysis for some values of the parameters q,p. We obtain 
r)v = 0.02294(20) at A 6 = and A 4 = 1.1, w = 0.02287(20) at A 6 = 1 and A 4 = 1.9, and 
rju = 0.02305(20) for spin-1 and D = 0.641. Taking again into account the uncertainty of 



A4 and D* we then obtain the estimate reported in Table LI I 



The CPRM was also employed in order to estimate the exponent a. It was applied to 
the 18th order series of x£, 2 /P an d Qa,o/P- The results are displayed in Table |XXIII| . We 
find ov = 0.0134(8) for A 6 = and A 4 = 1.1, ov = 0.0134(9) for A 6 = 1 and A 4 = 1.9, and 
ov = 0.0127(6) for spin-1 at D = 0.641. 



4. Ratios of amplitudes 

In the following we describe the analysis we employed in order to evaluate universal ratios 
of amplitudes, such as g 4 , r 2 j and q, from the corresponding HT series. In the case of g, c 2 , 
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TABLE XXIII. Results 


for 


a obtained using 


the CPRM applied to m 2 //3 and g 4 ,o//3 (18 


orrlprs i l~Tprp wp nspH t™;^ = 


0.9 


y 1 1 T* 

max x,x 5 J, max 


= 0.1. r-app is explained 


in the text. 






approx 


r-app 




A 6 =0, A4 = 1.08 




nTA2^> i r, n 

u-i-j *-±jq — 1 ,p — u 


29/34 


0.0134(8) 


A 6 = 0, A 4 = 1.10 




h A z^. 1 ^ n 


29/34 


0.0134(8) 






h A 1 m n 


53/62 


0.0133(10) 


A 6 = 0, A4 = 1.12 




hTA2« i n 

l'-lj i.ijg — -L)P — U 


28/34 


0.0135(9) 


A 6 = 1, A 4 = 1.86 




bIA2q = i >p= o 


29/34 


0.0133(9) 


A 6 = 1, A 4 = 1.90 




bIA2 9= i iP= o 


29/34 


0.0134(9) 






bIA2,j=2,p=o 


52/62 


0.0132(12) 


A 6 = 1, A 4 = 1.94 




bIA2 ?= i lP= o 


28/34 


0.0135(9) 


spin-1, D = 0.633 




bIA2 9= i iP= o 


21/34 


0.0127(5) 


spin-1, D = 0.641 




blA2 g =i, p =o 


21/34 


0.0127(5) 






bIA2 q= 2, p =o 


37/62 


0.0128(6) 


spin-1, D = 0.649 




b!A2 9= i >p= o 


21/34 


0.0126(5) 



C3 and C4 we analyzed the series P 3 ^ 2 g4 = X^lo a */^; P~ 4 °2 = Si=o a */^; P~ 3 °3 — Y.Ho a iP\ 
and /3~ 2 c 4 = Si=o a $ % ■ m or der to obtain estimates of the universal critical limit of g 4 , r 2 j, 
and Cj, we evaluated the approximants of the corresponding HT series at (3 C (as determined 
from the analysis of the magnetic susceptibility), and multiplied by the appropriate power 
of i. 

For an rath order series we considered three sets of approximants: Pade approximants 
(PA's), Dlog-Pade approximants (DPA's) and first-order integral approximants (IAl's). 
(I) [l/m\ PA's with 

l + m>n-2, (A32) 
Max[n/2 - g,4] < l,m <n/2 + q, (A33) 

where I, m are the orders of the polynomials respectively in the numerator and denominator 
of the PA. The parameter q determines the degree of off-diagonality allowed. The best 
approximants should be those diagonal or quasi- diagonal. So we considered PA's selected 
using q = 3. As estimate from the PA's we take the average of the values at (3 C of the non- 
defective approximants using all the available terms of the series and satisfying the condition 
( |A33|) with q = 2. The error we quote is the standard deviation around the estimate of the 



results from all the non-defective approximants listed above. We considered defective PA's 
with spurious singularities in the rectangle defined in Eq. ( |A21|) with x min = 0.9 (x min = 
only in the case of r w ) a; max = 1.01 and y max = 0.1. 

(II) [l/m] DPA's with 

I + m > n — 2, 

Max [(n - l)/2 - g, 4] < I, m < (n - l)/2 + q, (A34) 

where I, m are the orders of the polynomials respectively in the numerator and denominator 
of the PA of the series of its logarithmic derivative. We again fixed q = 3. The estimate with 
the corresponding error is obtained as in the case of PA's. We considered defective DPA's 
with spurious singularities in the rectangle with x min = 0, x max = 1.01 and y max = 0.1, cf. 
Eq. 

(III) [mi/mo/k] IAl's given by Eq. ( |A16| ). The off-diagonality parameter was fixed to 
be q = 3, and p = 1. As estimate we take the average of the values at (3 C of all non- 
defective approximants listed above. The error we quote is the standard deviation around 
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TABLE XXIV. Results of PA, DPA, and IA1 analyses of the series for g±, r2j and Cj. When 
results are not reported, it means that for that quantity no acceptable results were obtained from 
that class of approximants. The fraction at subscript is the number of non-defective approximants 
over the total number of approximants. The last column contains the estimates obtained by 



combining the three classes of approximants. 









PA 


DPA 


IA1 


combined 


9*4 


A 6 = , A 4 

A6 = 1 , A4 
spin-1, D = 


= 1.10 
= 1.90 
= 0.641 


23.500(60) 17/ 2i 
23.487(45) 17/21 
23.486(19) 20/21 


23.491(25)i 6 /i 8 
23.474(46) 17/18 
23.492(88) 17/18 


23.504(18)49/73 
23.498(24) 57/ 73 
23.491(52) 53/73 


23.499(16+20) 
23.491(21+40) 
23.487(18+20) 


re 


A 6 = , A 4 
A6 = 1 , A4 
spin-1, D = 


= 1.10 
= 1.90 
= 0.641 


2.051(13) 19/21 
2.052(12) 18/21 
2.0493(65) 20/21 


2.058(12)n/i8 
2.063(14) u /i8 
2.0461(24)i 6 /i 8 


2.048(7) 32/73 
2.048(4)33/73 
2.0456(16)23/73 


2.051(7+2) 
2.050(5+4) 
2.046(2+3) 


rs 


A 6 = , A 4 
A6 = 1 , A4 
spin-1, D = 


= 1.10 
= 1.90 
= 0.641 


2-24(9) 17/18 

2.23(11) 18/18 

2.40(8) 16/18 


2.21(13)i 7 /2i 

2.23(9)i 7/ 2i 

2.31(5)i7/2i 


2.23(5)37 /69 
2.23(5)36/69 
2.42(13)2 6 / 69 


2.23(5+4) 
2.23(5+6) 
2.34(5+3) 


no 


A 6 = , A 4 
A6 = 1 , A4 
spin-1, D = 


= 1.10 
= 1.90 
= 0.641 


-14(5) 16/21 
-14(5) 14/21 
-10(21) 13/21 




-13.3(1.3) 6/ 6i 
-12(4)io/ei 
4(36)i4/6i 


-14(4+0) 
-13(5+0) 
-8(25+0) 


10 4 C2 


A 6 = , A 4 
A6 = 1 , A 4 
spin-1, D = 


= 1.10 
= 1.90 
= 0.641 


-3.582(8) 15/15 

-3.574(7)i4/i6 
-3.570(12)i 5/ i 5 


-3.580(29)i2/i2 
-3.574(26)i2/i2 
-3.562(36)u/i2 


-3.586(24) 24 /33 

-3.585(38) 24/ 33 
-3.554(28) 25/ 33 


-3.582(7+6) 

-3.574(7+20) 

-3.568(11+4) 


10 4 c 3 


A 6 = , A 4 
Ag = 1 , A 4 
spin-1, D = 


= 1.10 
= 1.90 
= 0.641 


0.087(8)i2/i5 
0.086(5)ii/i5 
0.095(14)i 4 /i5 


0.080(10) 3/ i2 
0.078(12) 2/ i2 
0.100(12) 2/ i2 


0.084(7)26/36 
0.086(5)26/36 
0.090(4)3 /3 6 


0.085(6+0) 
0.086(4+0) 
0.090(4+0) 



the estimate. We considered defective IAl's with spurious singularities in the rectangle and 
^min = 0, x max = 1.001 and y max = 0.1. 

As in the case of the critical exponents, sometimes we also eliminated seemingly good 
approximants whose results were very far from the average of the other approximants. In or- 
der to arrive at a final estimate, the results from PA's, DPA's, and IA's were then combined 
taking also into account the relative number of non-defective approximants (before combin- 
ing the results we divided the apparent error of each set of approximants by the square root 
of the ratio between the number of non-defective and the total number of approximants). 
Of course all the above procedure to arrive at a final estimate is rather subjective. But we 
believe it provides reasonable estimates of the quantity at hand and its uncertainty. We 
report in Table |X.X1V| the results of each set of approximants, so that the readers can judge 
the reliability of our final estimates. The second error in the combined estimate is related 
to the uncertainty of the value of A4 and D*; it is estimated by varying A4 in the range 
1.08-1.12 for A 6 = 0, 1.86-1.94 for A 6 = 1, and D in the range 0.633-0.649 for the spin-1 
model. Errors due to the uncertainty of (3 C are negligible. 

We have also performed analyses of the series of for A6 = and several values of 



A4 by employing the Roskies transform JTTJ]. The idea of the Roskies transform (RT) is to 
perform biased analyses which take into account the leading confluent singularity. For the 
Ising model, where A ~ 1/2, one replaces the variable /3 in the original expansion with a 
new variable z, defined by 1 — z = (1 — P/fic) 1 / 2 . If the original series has square-root scaling 
correction terms, the transformed series has analytic correction terms, which can be handled 
by standard PA's or DPA's. Note that in principle IAl's should be able to detect the first 
non-analytic correction to scaling, but they probably need many more terms of the series, 
and practically need to be explicitly biased as in the case of PA's and DPA's. Indeed the 
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TABLE XXV. Details of the analysis of the 17th order series for /3 — 3//2 (74 with and without 
the use of the RT for some values of A4 and A6 = 0. In the PA, DPA and IA analyses with RT we 
used q = 2 (other approximants turned out to be much less stable). We fixed x m \ n = 0, x max = 1.1 
and y m ax = 0.25 for PA and DPA, and cc m i n = 0, x max = 1.01 and y max = 0.25 for IA. In order 
to perform a homogeneous comparison we used the same procedure for the direct analysis without 
RT (except that we used x max = 1.01 and 

£/max — 0.1). The fraction at subscript is the number of 
non-defective approximants over the total number of approximants. 



A 4 




PA 


DPA 


IA 


combined 


0.5 


direct 


22.62(58) 10/15 


22.43(16) 9/12 


22.55(19) 19/37 


22.48(15) 




RT 


23.75(27) 8/15 


23.48(25) 8/12 


23.29(50) 15/37 


23.56(23) 


0.7 


direct 


23.04(28) 10/15 


22.92(14) 9/12 


22.95(13) 15/37 


22.94(12) 




RT 


23.62(31) 13/15 


23.54(26) 9/12 


23.35(35) 17/37 


23.54(20) 


1.0 


direct 


23.58(23) 9/15 


23.40(6) 8/12 


23.38(19) 15/37 


23.41(7) 




RT 


23.56(27) 13/15 


23.57(15) 9/12 


23.45(21) 14/37 


23.55(14) 


1.1 


direct 


23.500(64) 12/15 


23.491(23) 11/12 


23.494(16) 17/37 


23.493(16) 




RT 


23.56(22) 13/15 


23.59(16) 9/12 


23.44(17) 14/37 


23.55(13) 


1.2 


direct 


23.63(4) 10/15 


23.613(9) 11/12 


23.612(9) 20/37 


23.613(8) 




RT 


23.56(20) 13/15 


23.54(34) 11/12 


23.43(16) 12/37 


23.52(15) 


1.5 


direct 


23.93(4) 14/15 


23.92(6) 7/12 


23.93(4) 




RT 


23.55(27) 12/15 


23.53(31) 11/12 


23.41(13)ii/ 37 


23.48(16) 


2.0 


direct 


24.14(28) 15/15 


24.07(17) 3/12 


24.15(9)31/37 


24.15(9) 




RT 


23.61(22) 12/15 


23.52(18) 10/ i 2 


23.44(12) 17/37 


23.50(12) 


3.0 


direct 


24.42(14) 15/15 


24.14(39) 6/12 


24.40(19) 15/37 


24.40(12) 




RT 


23.61(23) 14/15 


23.47(14) 11/12 


23.34(18) 16/37 


23.48(11) 


00 


direct 


24.78(10) 15/15 


24.57(19) 10/12 


24.81(16) 9/ 37 


24.75(9) 




RT 


23.59(20) 13/15 


23.47(11)iq/i2 


23.48(16) 13/37 


23.50(10) 



IA1 results without the RT turn out to be substantially equivalent to those obtained using 
PA's and DPA's. In Table XXV we report the details of the analysis without and with RT 
for some values of A4 and A6 = 0. These results are plotted in Fig. |I[ 



APPENDIX B: UNIVERSAL RATIOS OF AMPLITUDES 



1. Notations 



Universal ratios of amplitudes characterize the behavior in the critical domain of ther- 
modynamical quantities that do not depend on the normalizations of the external (e.g. 
magnetic) field, order parameter (e.g. magnetization) and temperature. Amplitude ratios of 
zero-momentum quantities can be derived from the critical equation of state. We consider 
several amplitudes derived from the singular behavior of: the specific heat 

C H = A ± \t\~ a , (Bl) 

the magnetic susceptibility 

X = C ± \t\~\ (B2) 
the spontaneous magnetization on the coexistence curve 

M = B\t\~ p , (B3) 
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the zero-momentum connected n-point correlation functions 

Xn = C±\t\-^ n ~ 2 ^ 8 . (B4) 
We complete our list of amplitudes by considering the second-moment correlation length 

e^rtr, (B5) 

and the true (on-shell) correlation length, describing the large distance behavior of the 
two-point function, 

fgap = /gapl*!""- ( B6 ) 

One can also define amplitudes along the critical isotherm, e.g. 

x = C c \H\-^ ps , (B7) 
£ = f c \H\-" /l3S , (B8) 

u = f c g J H r m - (B9) 

2. Universal ratios of amplitudes from the parametric representation 

In the following we report the expressions of the universal ratios of amplitudes in terms 
of the parametric representation (|7.5|) of the critical equation of state. 
The singular part of the free energy per unit volume can be written as 

T = h m R 2 ~ a g(9), (BIO) 

where g(B) is the solution of the first-order differential equation 

(1 - 9 2 )g'{9) + 2(2 - a)6g{6) = (1 - 9 2 + 2p9 2 )h{9) (Bll) 

that is regular at 9 = 1. One may also write 

x -^- R92{9) > 92{d) - — (i - v + W) — ' ( } 

X3 ~14 93(61 M ~~ g 2 mi-e 2 + 2(3e 2 ) ' (B13) 

v m °n-^a(f)) (l-9 2 )g>(9)-2(2 1 + f3)9g 3 (9) 

X ' = 14 R M ' Ud) - g 2 (9)(l-0 2 + 2m ' (B14) 

Using the above formulae, one can then calculate the universal ratios of amplitudes: 
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s- 7 02(0) 



c/ 2 = 

«4 = 

R- = 



A+ 

T- 

91 

c- 

ct_ 
ci 

_ aA+C + 
B 2 

aA-C- 



9 2 

° ' 92{V0) 

1 \ -37-2/3 04(0) 



_ j\ -3 T -2/3 

-a(l-a)(2-a)( ( 9 o 2 -l) 2 X- 2 0(O), 



u 3 = R 3 



B 2 U U 2 ' 
_ C^B 



(c- 



C7B 2 



v 4 = 



R A 



{C-f 



+ 3^ = 0o92(0of Sg 2 (9 )g 3 (9 ) 2 - ^o) 



CtB 2 



-6 



TT100 

^0 



R 



C+B 



6-1 



x (5C C Y 
lim z- 5 F(z) = p 1_5 /i(l). 



(B15) 

(B16) 

(B17) 

(B18) 
(B19) 

(B20) 

(B21) 

(B22) 

(B23) 
(B24) 



Using Eq. (7/7) one can compute F(z) and obtain the small- z expansion coefficients of the 
effective potential r 2 j in terms of the critical exponents and the coefficients h 2 i+i of the 
expansion of h{6). 



APPENDIX C: APPROXIMATION SCHEME FOR THE PARAMETRIC 
REPRESENTATION OF THE EQUATION OF STATE BASED ON 

STATION ARITY 

The parametric form of the critical equation of state, described by Eqs. ( |7.5| ), ( |7.6| ), and 
( [7.7|) , shows a formal dependence on the auxiliary parameter p. 

However all physically relevant amplitude ratios are independent of p, because they may 
be expressed in terms of the invariant function F(z) and its derivatives, evaluated at such 
special values of z as z = 0, z = oo and z = z , where F(z ) = 0. Notice that, despite 
the apparent dependence generated by the relation z = z(p,9 (p)), from the definition it 
follows that z must necessarily be independent of p. 

We can exploit these facts in order to set up an approximation procedure in which the 
function h(p, 9), entering the scaling equation of state, is truncated to some simpler (polyno- 
mial) function (p, 9) and the value of p is properly fixed to optimize the approximation. 

We found that, at any given order in the truncation, it is possible and convenient to 
choose p in such a way that all the (truncated) universal amplitude ratios be simultaneously 
stationary against infinitesimal variations of p itself. 

Starting from h^(p,9) we may reconstruct the function 
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(1 _^2)/3+ 7 - 



(CI) 



In order that all truncated amplitudes be simultaneously stationary in p, it is necessary that 
the function F^'(p, z) = F^\p, 9(p, z)) be stationary with respect to variations of p for any 
value of z. 

We shall prove that for any polynomial truncation h^(p, 9) it is possible to find a value 
p t , independent of z, such that 



dF®(p, 



dp 



0. 



(C2) 



p=pt 



a property which we shall term "global stationarity" . 

In order to prove our statement, let us rephrase the above condition into the form 

dF®(p,9) , dF®(p,9) d6_ n 

dp + 36 dp~^ {LS} 

where the implicit function theorem allows us to write 



de_ 

dp 



The definitions (715) and ( |7.7| ) imply 

dz z 
dp p 

Moreover it is trivial to show that 



dz/dp 
dz/dl 



dz (I 

de = z [e + 



(C4) 



2(39 
1-0* 



(C5) 



op 



'1 



dh®(p,ey 



p h®{p,B) dp 



dF®{p, 



dd 



i-e 2 



+ 



dh®{ P ,ey 



h®(p,6) d6 



(C6) 
(C7) 



Substitution of these expressions into Eq. ( |C3| ) leads to the following form of the global 
stationarity condition: 

[1 - (1 + 2 7 )^«(p, 9) + [1 + (2/3 - mp dhit) d ^ 9) - [1 - fl ^y* = 0. (C8) 

Let us now write down h®(p, 9) as a power series in the odd powers of 0: 

t-i 



h( t \p,9) = 9 + Y,h2n + i(p)9 2n+1 . 



(C9) 



n=l 



The series-expanded stationarity condition then takes the form 
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t-l 

E 

n=l 



9 



h 2n +i{p) + 



(2/? - l)p^ - 2 7 + 2n - 2 



h 2n -i(p) \ 9 



2n+l 



(CIO) 



with the convention h± = l. 

Let us now notice that, in the absence of truncations, the above equation must be 
identically true, since the original function F(z) is totally independent of p. This fact 
implies that the coefficients of the above power-series expansion must vanish individually, 
and this gives us an infinite set of recursive differential equations for the functions h 2n+ i(p), 
which must be automatically satisfied when the coefficients h 2n+ i(p) are properly defined. 

A truncation corresponds to arbitrarily suppressing all coefficients starting from h 2t+ i(p). 
Hence the global stationarity condition simply amounts to requiring 



(2/3 - l)pj- - 2 7 + 2t - 2 



h 2t -\{p) 



l2t+l 



0. 



(Cll) 



because all other terms vanish. The resulting equation can be solved by choosing p t such 
that the term in curly brackets vanishes, independent of 9. This concludes our proof. 

The effectiveness of this scheme is beautifully illustrated by its lowest-order implementa- 
tion, corresponding to the so-called "linear parametric model" , in the context of the three- 
dimensional Ising model. 

Let us truncate the exact scaling function h(p, 9) to its cubic approximation 



h^(p,9) = 9 + h 3 (p)9 3 



(C12) 



where h 3 (p) is taken from Eq. ( |7.13| ). Substituting h$ into the stationarity condition ( |C11|) 
for t = 2 we find 



i(2/3- 7 )p 2 -2 7 (l-7) 



0. 



which leads to 



P2 



/6 7 ( 7 - 1) 



/ 7-2/3 

The truncated scaling function vanishes at the value 6q\ 



(C13) 



(C14) 



91 



1 



7-2/9 

h 3 (p 2 ) 7(1-2/3) 



(C15) 



In this approximation the scaling equation of state turns out to be expressible simply 
in terms of the critical exponents f3 and 7. As a consequence, all the universal ratios may 
then be approximated to lowest order by appropriate algebraic combinations of the critical 
exponents. 

The above results reproduce the old formulae by Schofield, Litster, and Ho p5] , who 
obtained expressions for critical amplitudes in terms of critical exponents from a minimum 
condition imposed on the predictions extracted from a parametric scaling equation of state. 
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In the case of the linear parametric model, the globality of the stationarity property 



introduced by the above authors was shown by Wallace and Zia |27j , who adopted a slightly 
different, but essentially equivalent, formulation of the above model. 

As we showed above, global stationarity can be imposed on parametric models regardless 
of the linearity constraint. The next truncation, corresponding to t — 3, can also be treated 
analytically. Our starting point will be 



h^\p,9)=9 + h 3 (p)9 3 + h 5 (p)9 5 . 



(C16) 



The coefficients h^ 3 \p) and h^\p) are reported in Eqs. ( |7.13| ) and ( |7.14| ). By applying the 
stationarity condition (|C11| ) to h Fj (p), we obtain 



P3 



(7 -2(3) (1-7 + 2/3) 
\ 12(4/3 - 7 )F 5 



1 - 



\ 



72(2- 7)7(7 -l)(4/?- 7 )F 5 
(7-2/3) 2 (l -7 + 2/3) 2 



(C17) 



The truncated scaling function vanishes when 9 takes the value 6*0, which is now given 
by the relation 



01 



h 3 {p 3 



4h 5 (p 3 



2h 5 (p 3 ) \\ H(p 3 ) 



(C18) 



A general feature of truncated parametric models is the possibility of making a prediction 
on higher-order coefficients i^n+ij for n > t, in terms of lower-order coefficients. This is a 
natural consequence of having included by the parametrization some information on the 
asymptotic behavior of F(z) for large z. In practice we observe that each F 2n +x appears 
first in the coefficient h 2n+ i, in the form of a free constant of integration in the solution of 
the recursive differential equation relating h 2n+ i to h 2n -i- Since a truncation corresponds to 
setting h 2n +i = starting from n = t, this fixes the (truncated) value of all F 2n+ i starting 
from F 2t+1 . 

As an important consequence of this mechanism we observe that truncated models devi- 
ate from the exact solution only proportionally to the difference between exact and predicted 
coefficients, and this difference may be quite small even for very low-order truncations. 

In order to turn the above considerations into quantitative estimates, we need to get 
further insight into the properties of the functions h 2n+ i(p), especially in the vicinity of the 
stationary point p t . To this purpose we introduce the expansion 



h 2n +i(p) = 



2m 771 



2m+l- 



(C19) 



m=0 



Substituting this expansion as an Ansatz into the recursive differential equations, we check 
that F 2m+ i act as free parameters (integration constants), while the coefficients c„ )m must 
obey the following algebraic recursive equations: 

(n - m)c n)m = {(2(3 - l)m - 7 + n - l]c n _i )TO , (C20) 

for all n > m, subject to the initial conditions c m ^ m = 1. It is possible to find a closed-form 
solution to Eq. ( |U2T?| ): 
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{n-mj. k=l 



I n—m 

H(2/3m-j + k-l), (C21) 



but for our purposes the recursive equations will sometimes be more useful than their explicit 
solutions. 

Let us define the coefficients of the z-expansion of the truncated scaling function evalu- 
ated at the stationary point: 



m=2 



By definition, i^m+i coincide with their exact value i^m+i f° r & U m < t, while for m > t 
they are determined by the condition h^m+iiPt) — 0, which, according to Eq. ( |C19|) , implies 



J] c n , mP 1 m F^ +1 = (C23) 

m=0 



for all n > t. 

We can now prove the following lemma: 



^mc n , mPt 2m F 2 ( 2 +1 = (C24) 

m=l 



holds for all n > t. 

The proof is by induction. Let us assume the lemma to hold for a given value n; then, 
as a consequence of Eqs. ( |C24| ) and ( |C23| ), we obtain 

n 

E [(2/? - l)m - 7 + n]c ntm p 2 t m F^ +1 = 0. (C25) 

m=0 

Notice that the above equation holds also for the initial value n = t — 1, since in that case 
it coincides with the global stationarity condition. 



By use of the recursion equations ( |C20| ) we now obtain 

n 

E (« + 1 - m)c n+ i, mP l m F^ +1 = 0. (C26) 



m=0 



Because of the factor (n + 1 — m) the sum can trivially be extended up to n + 1, hence: 



71+1 71+1 

(n + 1) E c n+1 , mP 2 t m F^ +1 = E c n+1 , mPt 2m F£ +1 . (C27) 

m=0 m=l 

The l.h.s. vanishes by definition (cf. Eq. ( |C23| ) ) , hence the r.h.s. vanishes and the proof is 
completed. 

The above lemma is instrumental in evaluating the difference between the predictions 
originated by two subsequent truncations. By applying once more the definition of -F 2 m+i 



one can easily show that 
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E C "-,m Pt+liHm+l ~ Hm+l) + (Pt+l ~ Pt^F^m+l ~ (C28) 



m=0 

for all n > t. Let us now expand the equation to first order in the difference pl +1 — p%, and 
make explicit use of the lemma to obtain 

n 

E C n,mPt m (-^2m+l _ F 2m+l) — (C29) 
m=t 

for all n > t. It is crucial that i^m+i ~~ ^im+i = f° r an m < t. 

The equation we obtained allows us to express (within the approximation) all differences 
i^im+i — -^im+i m terms of the single quantity 

8F t = F 2t+1 -F^ +1 . (C30) 

Knowledge of the c n ^ m and some ingenuity lead to the explicit solution of Eq. QC29|) : 

p (t+l) _ p (t) , $ Ft /poi\ 

r 2m+l r 2m+l — a t,m 2(m-t) \^ 01 ) 

Pt 

where, for all m > t, 

I i \m—t m—t—1 

dt, m = } _,^ im-l) II (2/3m-j-k), (C32) 
\ m v- k=i 

and obviously d tt t = 1- 

As a corollary of this result, by comparing Eq. (|C29| ) when t = 1 to Eq. (|C24j) when 

(2) 

t = 2, we may write down a closed-form expression for all i^m+i coefficients (m > 1): 

El (2) _ ^ fT"}*^ 
2m+1 ~ 6m^-i)' 1 j 

which completes our analysis of the linear parametric model. 

One may also show that SF t is related to the variation of p t by the (linearized) relation 

SF t = { E m^mpf ^F^) (pi, - p\). (C34) 



Wt ~ 7 V, 



m=0 



Our numerical estimates, presented in Table [XTI|, show that p\, x —p\ is indeed small (< 0.01). 



In order to evaluate amplitude ratios, as shown in App. ||, we must also reconstruct the 
functions g{9) and g2{9), by solving Eqs. (|B11| ) and ( |B12|) respectively. These functions may 



be expanded in even powers of 9, with coefficients that are functions of p satisfying the same 
differential equations as h,2 n +i, Eqs. (|C10|) . One may show that, for any given truncation 
(p, 9) and arbitrary values of p, 

g (t) (p, o) = E I £ c ™p 2m ^f~ 2 ° 2n+2 + A ^ - ° 2 ) w+ ^ ( C35 ) 



n=0 \m=0 
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where A is an integration constant reflecting the arbitrariness in the zero-field value of the 
free energy. One may also show that for n > t 

s^^~^)>a w+7 -* ) <c36) 

where the terms in the r.h.s. are the coefficients of the Taylor expansion for (1 — 9 2 ) 2/3+7 . 
As a consequence the constant A may always be chosen such that g^(p,9) is truncated to 
0(6 2t ) for any arbitrary choice of p. 

In turn one may also prove that, for any h^'(p, 9), 



n=0 



.m=0 



9 2n . (C37) 



Now, according to Eqs. ( |C23| ) and ( |C24| ) , when we choose for p the globally stationary value 
p t , the coefficients in square brakets vanish for all n > t. As a consequence for any t the 
value p t insures the truncation of gz\pt,8) to 0(9 2t ~ 2 ). Thus a unique feature of p t is the 
simultaneous and consistent truncation of h(9) and g2(9). 

Notice that we might start by imposing a global stationarity condition directly on a 
truncated g<i{.P, &), obtaining a different stationary value for p, and make use of Eq. ( |B12| ) in 
order to reconstruct the corresponding h(6). However in this case, since h(9) must be an odd 
function of 9, there is no arbitrary integration constant (which is physically a trivial conse- 
quence of the definition of a reduced temperature) and therefore h{9) cannot be truncated. 
The resulting parametric model is mathematically consistent, but in practice unappealing, 
because the calculation of 9q from the equation h(9o) = and the evaluation of universal 
amplitude ratios becomes quite cumbersome. 

The above described formalism can be usefully employed in the context of the e-expansion 
of the critical equation of state. Comparison with e-expansion results will also shed further 
light on the meaning and relevance of the results derived by the prescription of global 
stationarity. 

Our starting point will be the result of Wallace and Zia ||27|| , who showed that, when 
appropriate conditions are imposed on the zeroth order approximation, the parametric form 
of the critical equation of state is automatically truncated in the powers of 9 2 when expanded 
in the parameter e = 4 — d. For easier comparison, notice that the parameter b introduced 
by Schofield [E3] and used by Wallace and Zia is the same as our 9q, and the variable change 
from #0 to p poses no conceptual problem. 

In our reformulation one may state that, within the e-expansion, it is possible to choose 
to lowest order a value p in such a way that, expanding the parametric equation of state 
in 9 and e, one finds, for all n > 2, 

W(Po) = 0(e" +1 ), (C38) 

and this property should survive the replacement po —* Po + 0(e). 

As a first application of our formalism, we can verify the consistency of the above state- 
ments, by checking that, for all n > 2, the condition 
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E c n ,mPl m F 2m+1 = 0(e" +1 ) (C39) 



m=0 



implies 

Y,mc r , mP l m F 2m+l = 0{e n ). (C40) 

m=l 

The proof is by induction. Assuming the property to hold for a given n, and exploiting the 
fact that 2/9 — 1 = 0(e), we obtain 

E [(2/9 - l)m - 7 + n]c n , m p^F 2m+1 = 0(e" +1 ). (C41) 

m=0 

The initial condition, corresponding to the case n — 1, has the explicit form 

7(7-l) + ^(2 / 9-7)^ = 0(e 2 ), (C42) 

and is a definition of p . Notice that p = limp 2 , and in the Ising model p\ = 2. 
By applying the recursion equations we then obtain 

n 

E (« + 1 " ^)cn+i,mPo m ^+i = 0(e" +1 ). (C43) 

m=0 

The sum can trivially be extended to n + 1 and, recalling the hypothesis, we obtain 

n+l 

E mc n+ltm pl m F 2m+1 = 0(e n+1 ), (C44) 

m=l 

thus completing the proof. Along the same lines it is straightforward to prove that, for all 
n > 2, 

n 

E m k c n , m p 2 m F 2m+1 = 0{e n ' k+1 ) (C45) 

m=l 

for all integers k < n. The initial condition (n = k) is trivially satisfied for all n > 2: 

E m n c rhm p 2 m F 2m+1 = 0(e). (C46) 

m=l 

As a consequence the more general statement 

h 2n+1 (p) = 0(e« +1 ) (C47) 

holds for all p admitting an e-expansion and possessing the limit limp = po- 

This relation implies in turn that, expanding in e the coefficients F 2m+ i for m > 2 
according to 
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2m+l 



£/m*e fc , (C48) 



k=l 



when the f m k for m < k are known then all f m k for m > k are fully determined. 

As a simple application of the above we obtained the following closed form result: 



fmi = ) ^ n Po 2m lim — , (C49) 

where 7 = 1 + |e and p = \/2- 

Let us now consider the linear parametric model with global stationarity in the context 
of the e-expansion: p 2 satisfies the condition p 2 = po + 0(e), though it does not coincide 
(and is not expected to) with the e-expanded p value adopted by Guida and Zinn- Justin 



Now notice that for any higher-order truncation the stationarity condition is still solved 
by Pt = Po + 0(e), as shown explicitly by the above derived Eq. QC41|) . As a consequence, 



any stationary truncation is an accurate description of the e-expanded parametric equation 
of state up to 0(e t ) included. Actually the freedom of choosing p leaves such an expansion 
highly under-determined, and many other prescriptions might work, including that of fixing 
p (or alternatively 60) to its zeroth order value. It is however certainly pleasant to recognize 
that our approach based on stationarity falls naturally into the set of consistent truncations. 
As a side remark, notice that all the coefficients of the e-expansion of p t will in general be 
changed order by order in t, and will also in general be complex numbers. This fact will by 
no means affect the real character of the expanded physical amplitudes, and will not even 
prevent p and 9q from taking real values in the actual three-dimensional calculations. 
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